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PREFACE 


This book is an outcome of and closely follows the out- 
line of a series of lectures presented at the Harvard Engi- 
neering School during the season of 1930-1931. 

The increasing importance of noise and vibration in the 
engineering work of the General Electric Company has 
brought to the author’s attention many interesting prob- 
lems, and has indicated a need for a convenient reference 
work on this subject. The purpose of this book is, there- 
fore, to present in concise form a resume of the various 
aspects of vibration prevention in engineering, which have 
come directly within the author’s experience. 

The book is intended to serve as a reference work for 
practicing engineers and particularly for the students in 
the Advanced Course in Engineering of the General Elec- 
tric Company. 

No one can do effective work in this field without some 
grasp of the fundamentals on which the subject is based. 
The method of presentation has, therefore, been to develop 
each topic in a way which leaves no questions as to the 
basic logic involved, but at the same time to omit long 
analytical developments not essential to this understanding. 

Discrimination has been used in selecting and develop- 
ing only those aspects of elasticity and dynamics which 
in the author’s experience have been found to be particu- 
larly important in vibration work. 

Although many fields have been touched upon, which 
could have been much further expanded, the use of material 
and viewpoints already conveniently available in other 
text books has, so far as possible, been avoided. 

The scope of this volume comprises all the material 
necessary for handling effectively any usual vibration 
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problem. An cfTort has been made to enhance the value of 
the presentation by clear visualization of the theory- and by 
directness of attack on practical problems. 

The sections on the application of complex quantities 
to the solution of vibration problems represent methods 
developed and used by my associates E. H. Hull and 
\V. C. Stewart, to whom I wish to express m3' indebtedness. 

I also wish to express mj' appreciation of the assistance 
given me by Dr. B. L. Newkirk, and bj' Messrs. P. L. Alger, 
E. L. Thearle and A. L. Ruiz, who reviewed and criticized 
the manuscript. 

I desire particularly to acknowledge the interest and 
vision of Prof. Arthur E. Norton, of the Har\'ard Engineer- 
ing School, who was primaril}' responsible for the presenta- 
tion of the lectures of which this book is the outcome. 


ScilEXECTADY, XeW YoRK 
June, 1932 


A. L. Kimball. 
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CHAPTER I 

VIBRATORY MOTION IN GENERAL 

1. Vibratory Motion. — ^Any motion which repeats itself 
through approximately equal intervals of time is called a 
periodic motion. Periodic motion is a very common form 
of motion in nature, as a little consideration will show. 
The planets perform periodic motions around the sun. 
Indeed, the rotation of the earth on its axis is a form of 
periodic motion. The swing of the pendulum, the recipro- 
cating motion of a piston, or any motion classed under the 
head of vibrations is a periodic motion. 

A vibration may be defined as a periodic back-and-forth 
or reciprocating motion which reverses its direction twice 
every cycle. Vibrations fall into two general classes: the 
first, a free vibration which has no driving force but which 
when once started repeats itself without the aid of any 
external agency, usually dying out only gradually; and 
the second, a forced vibration which is maintained by a 
driving force of any frequency. 

2. Conditions Necessary for a Free Vibration. — In order 
that a given mass shall freely vibrate, a displacement from 
its rest position must be opposed by a so-called restoring 
force. When the body is displaced from its rest position, 
work is done against this force whose amount depends'^upon 
the amount of the displacement. When the body is set 
free, this potential energy is released by giving motion to 
the body so that when its rest position is reached the energy 
is all kinetic. If there are no energy losses during this 
motion the kinetic energy of motion, as the body passes 
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through its rest position, is equal to the potential energy' 
imparted by the original displacement. This carries the 
body onward to a new displacement on the other side of its 
rest position until again all the vibrational energy becomes 
potential. This process will continue indefinitely, the 
\'ibrational energy alternating between potential and 
kinetic. In ordinary’ cases, such a freely vibrating body 
will gradually come to rest because of friction of differ- 
ent kinds Avhich dissipates the vibrational energy^. Such a 
^•ibration, where the body is acted ujion by the restoring 
force alone, is called a free vibration. The frequency of a 
free vibration has one particular value called the natural 
frequency of that vibrating system. 

The phenomenon of vibration decay, that is, the process 
of gradual falling off of amplitude, is characteristic of a 
free vibration, such, for instance, as the dying out of the 
swing of a pendulum. 

3. Forced Vibrations. — Forced vibrations arc of quite 
different character from free vibrations, and should be 
clearly differentiated from them. 

A forced vibration is an oscillation of a body maintained 
by an c.\ternal periodic driving force. The natural fre- 
quency depends upon mass and the restoring force and has 
only one value, for a simple system, whereas a forced vibra- 
tion may have any frequency, depending on the frequency 
of the e.vternal stimulus including the natural frequency. 
In fact, a body with no restoring force whatsoever acting 
upon it will take up a forced vibration. The quiver of an 
unbalanced motor, always in step with its rotational fre- 
quency, is an example. 

The natural frequency depends upon the physical 
properties of the sj’stem itself, having nothing to do with an 
axternal stimulus. 

4. Resonance. — \^Ten the external stimulus which pro- 
duces the forced vibration is c.\actly in step with the natural 
frequency, that is, the frequency of free vibration of the 
body, the amplitude of the vibration may become very 
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large indeed, and the \dbration is said to be in resonance 
with the stimulus. In fact, in the ideal case where no 
friction is present, the amplitude of a resonant vibration 
rvill become infinite. 

In actual cases, friction is always present, so that the 
resonant amplitude builds up until the rate of energy dissi- 
pation in friction becomes just as large as the rate of energy 
supply through the stimulus. 

A resonant vibration may be thought of as a sustained 
free vibration. Although sustained by an external stimulus, 
it is not forced in the sense that a frequency other than the 
natural frequency is maintained. If the stimulus were 
removed, the vibration would continue at the same fre- 
quency "with slowly dying amplitude. 

An example of resonance of electrical oscillations famil- 
iar to all is found in the tuning of radio sets. The natural 
frequency of the electric circuit in the receiving set is made 
to correspond with the frequency of the oncoming electric 
waves, so that large resonant electrical amplitudes are 
built up. 

A powerful station in close proximity will produce forced 
vibrations of sufficient intensity to be audible without 
tuning being resorted to. 

An alternating electric current is an example of a forced 
electrical oscillation. The resonant condition is dangerous 
because of local high voltages, and is avoided. 

5. Vibrational Energy . — h body vibrating freely at a 
constant amplitude is a system of constant energy, as no 
energy is received from or given up to the outside. Thus, 
during the periodic cycle of transfer of energy between 
potential and kinetic the sum of these two must always be 
a constant. ■ In fact, this knowledge affords an important 
means of determining free vibration frequencies to be dis- 
cussed later. The potential energy at the rest positions of 
the body at the ends of its swings must just equal the 
kinetic energy at the instant of its maximum velocity as it 
passes through its position of zero displacement. 
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Such a vibrating system with no dissipation present is an 
example of a conservative syslcw of forces, all the energy 
being conserved at every instant within the system, regard- 
less of liow it moves. 

6. Elasticity and Vibrations. — By far the most common 
source of potential energs' in a free vibration is the potential 
energy’ of elastic deflection, and a consequent clastic restor- 
ing force. In this treatment, clastic restoring forces will be 
dealt with almost exclusi\'ely. 

Another common source of potential energy is that due 
to gravitation, the restoring force being gravitational. A 
gravity pendulum is an illustration in point. Many \’ibra- 
tions are produced by restoring forces, both clastic and 
gravitational in character, such as a gravity pendulum ^ 
whose frequency is increased by a spring at the top. In 
elastic vibrations, particularly when the frequency is high, 
the gravitational force is negligible and so is not considered. 

Centrifugal force may furnish a restoring action, as 
illustrated by the increase in stiffness of long turbine 
buckets, due to this cause, in a turbine wlicel revolving at a 
high speed. If the buckets are short and stiff, this effect 
is small and may be neglected in comparison with the elas- 
tic restoring force. 

Common examples of clastic restoring forces are the 
deflection of springs of different materials, the bending of 
bars, or the compression and extension of bars, the com- 
pression of a gas or liquid. In fact, anj’’ conceivable elastic 
distortion of anj^ elastic material may furnish the restoring 
force for some type of vibration. 

7. Vibrations and the Principle of Minimum Energy. — 
One of the broadest general principles in mechanics is the 
principle of miniinnin potential energy, whereby any system 
of bodies or particles, however they may be connected, 
tends to assume a shape such that its potential energ^'^ is a 
minimum. A ball rolls down to a position of minimum 
gravitational energy’. A compressed spring tends to 
e.\pand, until its potential energy’ is a minimum. Two 
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springs compressed against each other adjust their deflec- 
tions until their resultant potential energy is a minimum, 
etc., regardless of the complexity of the system, A vibra- 
tion is the result of the effort of an elastic system to seek its 
position of minimum energy. It acquires a velocity in this 
effort which carries it beyond the position of minimum 
energy, and it then seeks it again and then again, through 
many repeating cycles, if the friction is small. The posi- 
tion of minimum energy is evidently identical with the 
rest position as previously defined. Indeed, any elastic 
system in stable equilibrium must have the capacity to 
vibrate back and forth through its position of stability if 
the friction is small. The test of stability is minimum 
potential energy such that any small displace- 
ment of the system from this position causes 
the energy to rise. It is for this reason that 
the capacity to vibrate is such a universal 
characteristic of material bodies in nature. 

8. Degrees of Freedom. — The number of 
degrees of freedom of a system capable of 
vibration is measured by the number of inde- 
pendent free vibrations, each with its own 
characteristic frequency, which that system 

•' on Reed. 

may assume. 

The weightless reed of Fig. 1 carrying the mass m has 
(omitting possible torsional oscillations), two degrees of 
freedom, one in the plane of the paper and one of higher 
frequency where the mass vibrates in a plane perpendicular 
to the plane of the paper with an edgewise deflection of the 
reed. If the reed is broad and thin the latter frequency is 
not readily excited, and we have practically a system of one 
degree of freedom. 

The mass on the weightless helical spring of Fig. 2 also 
represents a system of one degree of freedom if the sideways 
motion is prevented by the pair of rods on which it slides. 
Without the rods it may assume a pendulum swing in 
either one of two directions perpendicular to each other 
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with a restoring force due to gravity, as well as a torsional 
oscillation about any one of three perpendicular axes, 
giving it six degrees of freedom in all. 

A cylindrical mass supported on three springs as shown 
in Fig. 3 also has six degrees of freedom, four of which 
arc easily excited, namely, a sideways swing, an up-and- 
down bounce, a torsional frequency about a vertical axis, 
and a sidewa 5 ’s shearing oscillation. These four inde- 
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Fig. 2. — One Degree of Freedom. Fig. 3. — Six Degrees of Freedom. 


pendent frequencies can be readily obseiA’^ed in an e.xperi- 
mental model. 

From a mathematical standpoint the number of degrees 
of freedom of a system is measured by the number of indc- 
pendent space coordinates by which its motion must be 
specified. 

A continuous elastic body like a long deflected rod may 
vibrate in any of an infinite number of frequencies, and thus 
represents a system of an infinite number of degrees of free- 
dom. This is characteristic of any continuous elastic body, 
assuming its molecular structure to be of infinitely small 
dimensions. From a practical standpoint, however, only 
the lower frequencies are easily excited. These cases will 
be subsequently discussed in detail. 



CHAPTER II 


SIMPLE LINEAR VIBRATION ; HOOKE’S LAW AND VIBRA- 
TIONS; ELASTIC CONSTANTS 

9. Simple Linear Vibration— Take the case of a weight 
W vibrating up and down on the end of a weightless spring. 
The frequency is given by the formula 


/ 


\ W 


( 1 ) 


Itt 

where k = elastic constant of the spring 
= force per unit deflection. 

If a downward deflection d is produced by the weight W, 
evidently 


W 

k = ~, so that (1) reduces to 
d 


f 


7-v yd 


l-K^d ■\/ d 


cycles per second 


where 


g = 386 inches per second squared, 
and d is measured in inches. 


. . . ( 2 ) 


w 




Thus the frequency of a system like that 
of Fig. 4 can be at once calculated when 
the gravity deflection d is known. Consider 
the case of a heavy motor on a compara- ^“Simple 

tively light structure. By measuring the tion. 

small downward deflection of this structure 
produced by the motor, the up-and-down natural fre- 
quency of the system is quickly estimated by formula (2). 
Cases often arise in engineering where this formula is 
applicable. 
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10 . Derivation of Formula for a Linear Vibration. — 
Formula ( 1 ) is easily derived by an application of Newton’s 
second law of motion, / = iiia where / = force and a = 
acceleration. If the mass ?;/ is gi\'en a displacement y, 
from its rest position, the restoring force = — ky from the 
definition of k. It is negative because the force is directed 
towards the origin in the negative 3' direction. Thus 

f = ma = m-r— = — ky 
or 

d~v 

m-j^ + ky ==0 ( 3 ) 

d V 

If V = 3’o and -7^ = 0 when / = 0 , 
at 


the solution of this difl'ercntial equation is 

3’ = 3'o cos oij . . 

where w, = and yo = amplitude of 3'. 


( 4 ) 


Since w. = 2-/. 


’ in 


The s3'mbol /c is used to denote the free vibration fre- 
quenc3' as distinct from any other frequenc}' as might arise 
in the case of a forced or sustained vibration. 

It is seen in Fig. 5 that the solution of equation ( 3 ) as 
given by ( 4 ) is represented graphically b}'’ the projection of a 
vector of length 3-0 on the horizontal diameter of the circle 
described b\' the tip of this vector as it revolves at the con- 
stant angular velocit3’’ w. The important point to note 
here is that this is simple harmonic motion, as simple har- 
monic motion is defined in just this way. 

11 . Hooke’s Law and Harmonic Motion. — In the 
derivation of equation ( 3 ) the restoring force was taken as 
~ky. This assumes that Hooke’s Law is obe3'ed, that stress 
is directly proportional to strain, or force is directly' propor- 
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tional to displacement. Thus, whenever the restoring 
force is directly proportional to displacement, the resulting 
vibration is simple harmonic, or just harmonic as it is 
usually called. Since Hooke’s Law that stress is directly 
proportional to strain is to a close approximation univer- 
sally true of practically all the elastic materials of nature, 
vibrations involving elastic deflections are also universally 
harmonic in character.* This is fortunate, as a harmonic 
vibration is the simplest possible type of periodic motion 
from several standpoints. It is easily handled mathe- 
matically. Its frequency remains the same regardless of 
amplitude. In fact, the vibrations of all musical instru- 
ments must be closely governed by 
Hooke’s Law; other\vise the indi- 
vidual notes would change with the 
intensity of the vibration, which 
common observation tells us is not 
true. 

12. Superposition of Vibrations ; 

Law of Similar Systems. — A uni- 
versal law of harmonic vibrations 
is that for systems of two or more 
degrees of freedom, including con- 
tinuous elastic media with an in- 
definite number of degrees of freedom, any number of the 
characteristic free vibrations may be superimposed upon 
each other without affecting each other either as to their 
individual frequencies or amplitudes. For any one char- 
acteristic free vibration of a system of several degrees of 
freedom, all component parts vibrate exactly together, that 
is, they all have their maximum velocities at the same 
instant and all come to rest at the same instant. The 
phase of some parts may be reversed with respect to' other 
parts, however, but this is always an exact 180° reversal, so 
that all parts vibrate in step. The portions of the system 
which are reversed in phase are separated from those which 

* For a rigorous mathematical discussion of this point see reference J . 
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arc unrcverscd by stationary' points called nodes. The 
irregular appearance of the vibration of such a system is due 
to two or more of these characteristic vibrations taking place 
on top of each other, independently. Take as an illustra- 
tion the simple case of a system of two degrees of freedom 
shown in Fig. 6 (considering up-and-down vibrations only). 

The two modes of vibration of this system are easily 
excited independently. For the lowest or gravest fre- 
quency, all parts oscillate up and down together, with the 
maximum amplitude at the bottom and a node at the top 
at the point of attachment. For the second or higher 
mode of vibration all parts vibrate in step as before, but 
7 ni and 1112 vibrate towards and away from 
each other in opposite phase with a node in 
the spring between them, and also a node 
at the top as before. All parts of the system 
between these nodes vibrate in opposite phase 
to those below the lower node. If the system 
is given a downward impulse, both vibrations 
take place at once, which, with keen observa- 
tion, can be separated by the eye. 

The possibility of superimposing vibrations 

* VI. X A * 

Degrees of Free- IB this Way follows from Hooke's Law. The 
dom. stresses produced by one system of displace- 
ments add directly to those produced by a 
second system of displacements, each system producing 
its own stress components regardless of the presence of the 
other. In other words, the stress components produced 
by one mode of vibration at any instant are the same as if 
the other mode of vibration were absent. The same rea- 
soning can be extended to systems having any number of 
modes of vibration. 

A very' interesting law, called the law of similar systems, 
can be shown to be true, on the basis of the laws of me- 
chanics and principles of elasticity' just given. This may 
be stated as follows; TT 7 /eK hvo vibrating systcjiis arc made 
of the same material and arc exactly shnilar in dimensions. 
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though not of the same size, their fregtiencies of vibration are 
inversely proportional to their linear dimensions. 

13 . Harmonic Analysis. — It was first shown by the 
French mathematician Fourier that any wave shape what- 
soever can be exactly represented by the so-called harmonic 
series 

y = Jo -{■ ji cos {oit 6i) 4 - y2 cos {2o}t -j- 62) 

4 - yz cos ( 3 w/ 4 " ^3) etc ( 5 ) 



fix) 

T" 

) 







It is necessary only to chose the amplitudes yi, j2, etc., 
and the phase angles 61, 62, etc., correctly. The individual 
terms in this series are called harmonics, because each by 
itself represents a simple harmonic motion. The lowest 
frequency term, containing co, is called the first harmonic 
or fundamental; the 
second one, containing 
la and of t%vice the 
frequency of the funda- 
mental, is called the 
second harmonic; the 
next the third harmonic, 
etc. The simple system 

of Fig. 4 has no higher harmonics, only the fundamental 
being present. The motion of any point of the S5'stem of 
Fig. 6 contains two frequencies. 

The periodic curves as shown in Fig. 7 are readily repre- 
sented by a series of the type ( 5 ), an infinite number of 
harmonics being required to obtain the mathematically 
exact curve. The series for the first of these is as follows: 

J W = -j sm h , sm h - sm h . . . 

ttI c 0 c 5 c 


Fig. 7. — Periodic Curves. 


If a few terms of this series are plotted against x as 
each successive curve thus obtained is added, the rec- 
tangular shape of the wave wll be seen to be more and 
more nearly approximated. 
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If a piano wire is struck a sharp sudden blow with the 
hammer, a wave with a relativelj'- sharp peak is started and 
travels back and forth in both directions. This wave 
contains a verj' large number of harmonics, the upper of 
which, however, soon die out, leaving the fundamental and 
lower harmonics as predominent, which give the tone its 
quality. 

14. Young’s Modulus, Stress and Strain. — Doubtless 
the first of the elastic constants to be defined was Young’s 
modulus, usually designated by the letter E. The reason 
for this is that it is comparatively easily evaluated from 
such c.\perimcnts as the force required to produce a certain 
elongation in a given wire or bar, being defined as the 
force per unit area of cross-section divided by the elongation 
per unit length produced in the stressed specimen. 

If force per unit area — F/ A — a, called the stress, 
and if the elongation of the bar = u = «/, where e = elonga- 
tion per unit length, called the strain, and I = length. 


F _u 
A ' I 



( 6 ) 


In other words, E is numerically equal to the stress neces- 
sarj’' to make u = /, or that required to double the length 
of the specimen, assuming that Hooke’s Law that stress 
is proportional to strain, is followed up to that point. 

The relation (6) which may be put in the form Et = a 
will be often referred to from now on, where a is called the 
tensile or compressional stress and e the strain. 

The stress <r may be thought of as the load on the faces 
of a unit cube of the stressed materia! perpendicular to 
the direction of the force, and e as the amount that the 
unit cube is elongated or compressed. 

IS. Shear Modulus. — If, besides Young’s modulus, 
the shear or rigidity modulus, designated by N, is known, 
the clastic properties of an isotropic material are com- 
pletely specified. An isotropic material is one whose phys- 
ical properties are the same in every direction. 
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Not only does a bar of material have the power to resist 
tension and compression, but it can also withstand torsion. 
Torsion sets up a stress which tends to produce a slipping 
or shearing of successive sections of the bar across each 
other. Shearing stress is defined as the stress per unit 
area which acts parallel to a surface instead of normal to 
it, and which tends to produce a slipping of successive 
layers. This causes a distortion of the material called a 
shearing strain. 

Figure 8 shows a unit cube whose faces perpendicular 
to the paper are acted upon by a shearing stress. Note 
that if one pair of opposite faces are so stressed, the other 
pair must be stressed in shear by the same ^ 

amount in order to maintain equilibrium. It 
is a general law of elasticity that perpendicular 
planes at any point in a material under stress 
must always have such equal balancing 
shearing stress components. If the shear on 
one plane is zero in all directions, it is zero in 
all perpendicular planes; or if zero in one direction it is 
zero in the corresponding perpendicular plane. 

Figure 8 shows the type of strain produced by a shearing 
stress on a unit cube. Its faces become slightly diamond 
shaped. 

The shear modulus N is defined in a way exactly 
analogous to Young’s modulus as given by equation (6), 
that is 




Fig. 8. — Shear 
on Unit Cube. 


N = 


( 7 ) 


where t is the shearing stress or load per unit area, and 
e, is the amount that the upper surface of the unit cube 
has slipped or sheared with respect to the lower surface. 
Note that e, is also equal to the angle of tilt of the vertical 
faces of Fig. 8 from the perpendicular. 

16. Poisson’s Ratio.— If a bar of elastic material is 
stretched or compressed by a small amount its diameter 
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docs not remam constant, but decreases slightly with ten- 
sion and increases with compression. The ratio of this 
sidewaj's contraction or e.xpansion of a unit cube to the 
corresponding stretch or compression is a constant for the 
material; it is called Poisson’s ratio, and it will be signified 
bj' the Greek letter Poisson’s ratio has a value wliich 
\-aries from I to J for ordinary elastic materials. It is never 
greater than I, otherwise the volume of the material would 
decrease with tension and increase with compression, which 
is impossible for an isotropic material. If clastic deforma- 
tion is attended bj' no volume change, Poisson’s ratio has 
the value of almost exactly T This condition is very'’ nearly' 
fulfilled by rubber and gelatin. Although E and N for 
these materials are small, they' are very' resistant to volume 
compression. An example of a material with Poisson’s 
ratio almost equal to zero is cork. A cylinder of cork can 
be compressed along its axis by' a large amount and show 
practically no change of diameter. Its elastic properties 
are rather poor, however, in that it does not immediately 
recover its original shape when the load is removed. 

The following table gives approximate values of Pois- 
son’s ratio for a number of materials: 


Table I 


Material 

Poisson's Ratio 

Iron 

0.28 

Steel 

0.28-0.31 

Copper 

0.35 

Brass 

0.33 

Glass 

0.24 

Rubber 

O.SO 

Cork 

0.00 


17. Relation between Elastic Constants. — ^The three 
elastic constants just defined are Young’s modulus E, the 
rigidity’ modulus jY, and Poisson’s ration v. Any two of 
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these constants completely determine the elastic prop- 
erties of a material, so there must be a definite relation 
connecting them. This relation may be found as follows: 

Let ABCD be the side of a unit cube before it is 
strained. Assume a tension a to be applied at the top and 
bottom faces so that the original shape ABCD of the side 
shown in Fig. 9 is distorted to the shape A'B'C'D. The 
tensile strain produced by o-is^^' = e. Poisson’s ratio 
is evidently C'C -f- A A' = v. It is seen from the figure 
that, although a produces no tangential stress on the top 
plane AB, it does produce a tangential or shearing stress 
across any plane inclined to AB. This shearing stress is a 
maximum on a plane tipped 45° from the horizontal, such 
as AC. It is understood that .4C is the edge of a plane 
perpendicular to the plane of the paper. The component 
of 0 - which is parallel to .4 C is evidently a -f- V^. Since 
plane AC has an area V2 times the unit area AB, the 
sheBx per unit area on AC produced by cr is (r/2, which is 
called T (§ 15). 

Therefore 

Shear on AC = r = ^, . (S) 

For clearness, the strains are 
shown greatly exaggerated in 
Fig. 9. Actually, they are so 
small that the change of areas 
and consequent change of stress 
intensities thus produced may be 
neglected. Such corrections are 
called second-order corrections 
and are universally neglected in 
all the mathematical theory of 
elasticity. It is only when the 
strains are so large that the shape of the body changes 
materially, as for a material like rubber, that this approxi- 
mation results in appreciable errors. Thus, in estimating 
the shearing strain produced by the stress r across plane 


Stress o cr 
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AC, the final length is taken equal to the original 
length A C, etc. 

( € V€ \ /“ 6 “1“ r6 

— 2 ~' 

But the shearing strain e, with respect to AC as defined in 
§ 15 is equal to 24>, so that 

2^ = € + J'€ = €, (9) 

where e = tensile strain and v = Poisson’s ratio. 

Using (6), (7), and (8) 

r = 5 = Y = Ne. or Ee = 2Ne. . (10) 

Substituting (9) in (10), Ee = 2Ne (1 + f), or finally 

E = 2N{\ + v) (11) 

This is the required relation connecting the three elastic 
constants of a material. 

Although here proved for one special case it can be 
shown to be generally true. (See reference 2.) 

18. Further Remarks on Elasticity. — It is seen that a 
single tension or compression always produces shearing 
stress in a material. The only possible type of stress in 
which shear is absent is found when the compression or 
tensile stress is the same in ever>' direction, as is the case 
at any point in a fluid under pressure. It is evident that 
a fluid cannnot resist a shearing force. The slightest 
shear on a fluid causes flow, due to the slipping over 
each other of adjacent layers. Thus, the shear modulus 
of a fluid is zero. Young’s modulus also has no meaning, 
in the case of a fluid, since, as just e.\plained, a stress in one 
direction must cause shearing forces which in turn cause flow. 
Thus the only elastic modulus possible for a fluid is the 
compressibility modulus or bulk modulus which is a 
measure of its resistance to hydrostatic pressure. The 
bulk modulus of air is small, but that of water is compara- 
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lively high. Of course, a solid also has a bulk modulus, 
but this is readily expressed in terms of any two of the 
three elastic constants previously defined. 

For further explanations of these quantities and how 
they are derived, works on elasticity should be consulted. 
(See reference 2.) The purpose of the discussion of elas- 
ticity just given is to furnish a basis for the discussion of 
vibrations which follows. Nearly all vibrations involve 
elastic deformations, so that an appreciation of some of 
the basic laws of elasticity as given is essential to the study 
of vibrations. 



CHAPTER III 


VIBRATION OF LOADED CANTILEVER; BENDING 
MOMENT AND POTENTIAL ENERGY AT SECTION 
OF DEFLECTED BAR 

19. Vibration of Weightless Cantilever and Beam with 
Load. — If a weight is placed on the end of a cantilever 
as shown in Fig. 10, the weight being several times more 
massive than the cantilever, the system may be treated as 
one of a single degree of freedom and the frequency found 
to a first approximation by neglecting the weight of the 
cantilever. 

The simplest way to estimate the frequency is to 

measure the gravity de- 

flection d in inches pro- 

— - ■ . ■ - I duced by the weight, 

^ measuring this in a line 

Fig. lo-Cantilcvcr with Load. Opposite the center of 

gravity of the weight. 
The frequency is then readily calculated by formula (2), 
that is, / = 3A3/\/d cycles per second. 

If d is not readily measured, the formula for d as pro- 
duced by a weight IF placed on the end of a cantilever 
may be used. This is d = Wl^/3EI, where I = length 
of the cantilever to the center of graA'ity of the weight, 
E is Young’s modulus for the cantilever material, and I 
the moment of inertia of its area of cross-section, some- 
times called the second moment of this area. If the cross- 
section is rectangular, I = hh^/\2, where b = breadth and 
h = thickness or depth. The value of d thus calculated 
may then be used in place of the measured value. 

If the elastic member is a beam supported by pivots 
at both ends, which carries the weight at its middle point, 

18 
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instead of a cantilever, the process is just the same except 
that for such a beam d = WP/^SEI. The method can 
evidently be extended to beams supported in other ways 
by using the right expression for d. (d = deflection at 
center of gravity of the weight.) 

It is evident that neglecting the weight of the beam 
leaves out of account some of the inertia of the vibrating 
system, so that the frequencies so calculated come out a 
few per cent too high. 

20. Elastic Constant of Cantilever and Beam. — The 
deflection d at the end of a cantilever is found from the 
general formula for bending moment in a bar 





Fig. 11. — Uniform Cantilever. 


where M = bending moment, I = the second moment of 
the area of cross-section, and d~yjdx^ is the curvature of 
the bar due to elastic deflection at the section under con- 
sideration. The complete derivation of this important 
formula will be given in a later section. Consider the 
cantilever of uniform cross-section shown in Fig. 11 car- 
rying a weight W on its end. Let y be the downward 
deflection from the horizontal at any point x. The bend- 
ing moment at x due to W = W{1 — oc). 

Thus 


Integrating: 


M = Ei-g = W(l - X) 


( 13 ) 


When dy/ dx = 0, x = 0, so that Ci = 0. 
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Integrating again: 



When y = 0, .x* = 0, so that C 2 = 0. Thus: 



This is seen to check the expression for d given in § 19. 
The elastic constant 




If this is substituted in (1) 


IF 3EI 


P 


f = ±JMi . 

■' 27r^ lF/3 


3.13 


/3l7 

\ w 


WP 


(16) 


(17) 


The same result would of course be obtained by substitut- 
ing (15) in (2) as indicated in § 19. 

The calculation of the elastic constant for a beam of 
uniform cross-section supported and loaded in other ways 
than that given can be calculated by the method just given. 

If the cross-section is not uniform, I is not a constant 
but depends upon x, which must be taken account of. 
The problem is of course more complicated in such a case. 

21. Derivation of Bending Moment Formula. — The 
formula 




just used is derived as follows. 

Figure 12 shows a cross-section of a bar subject to 
bending. Let the .v axis coincide with the neutral axis, 
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and the upper half be in tension and the lower in com- 
pression. 

It will first be shown that the neutral axis must pass 
through the center of gravity of this area. 

For equilibrium the total tension where y is positive 
must equal the compression where y is negative. 

Expressed mathematically, h n 


If ffdxdy = 


Unit length 


where a is the stress intensity, and 
the integration covers the entire 
area of cross-section. 

Consider a section of the rod of 
unit length, bent from the straight 
to a curve of radius of curvature 
as shown in Fig. 13. In actual 
cases, p is very great compared 
with the diameter of the rod. 


Tension 



7 = Radius 
of 

curvature 


Compression 


Fig. 12. — Cross-Section of Bar. 


Fig. 13. — Deflected Section 
of Rod. 


^ Consider a fiber EF at a distance y from the neutral 
axis as shown in Fig. 13. The elongation of this fiber is 
equal to the strain e, since the section is of unit length. 

If = the angle made by AB and CD extended. 


e = ydi = - 
P 
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The angle 0\, which is seen to be the angle of cur\-ature 
per unit length, is called simply the curv'^ature; its reciprocal 
is equal to p, the radius of this cur\'ature. In an actual 
unit length, p maj’’ change somewhat, as the bending 
moment varies along the unit length. Here we have 
assumed a slice of infinitesimal thickness to be expanded to 
unit length, all of which carries the same bending moment 
and thus has the same radius of cur\’’ature p. 

From equation (6), 


= Ee 


From equation (19) 


Therefore 


p£ = 


El 

E 


Ev 


a — 


( 20 ) 


Substituting (20) in (18), 

—J' J' ydxdy = 0 . . . 


. . ( 21 ) 


from which it is seen that the .v axis, taken as the neutral 
axis, passes through the center of gravity or centroid of the 
area of cross-section. By definition, a line must pass 
through the center of gravity of an area when the sum of 
ever)" element multiplied by its perpendicular distance to 
that line on one side of that line equals the sum of the 
same quantities taken on the other side of the line, which 
condition is satisfied by (21). 

From Fig. 12 the total moment about the neutral axis 
is seen to be 


M = J' J' aydxdy — J' J" dxdy = — . (22) 

where / = J' J" y~ dxdy integrated over the whole section. 


J is recognized to be the second moment of cross-section 
of the rod (§19), It is a geometrical property of the 
cross-section of the rod and has no physical correspondence 
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to the moment of inertia of a body usually specified by the 
same symbol. 

From calculus it is shown that 


1 

P 


d?y 



which equals 


dx^ 


when dyfdx is small as it is in deflecting beams when x is 
measured along the beam. Substituting this in (22) : 


M 



(23) 


which is the bending moment formula. 

22. Potential Energy of Bending. — If W i is the elastic 
potential energy of bending of a rod per unit length, 


Wi 


1 


In general, 


Work = W 


tOzY 

Vdx-V 

= ff,s 


where / = force and 5 = distance measured in the direc- 


tion of the force. 


For work of bending, 



where L = torque = bending moment M, and 6 = angle 
through which the rod is bent, measured in the direction 
of the torque. 

Thus in bending unit length of the rod through the 
angle d, 




El Odd = 


Eldi^ 





(24) 


The substitutions used are evident by referring to the 
previous section. 

This is a useful formula in certain types of bending 
vibrations. 



CHAPTER IV 


SIMPLE TORSIONAL VIBRATION; MOMENT OF INERTIA; 
VIBRATION WITH ONE DEGREE OF FREEDOM AND 

DAMPING 


23. Simple Torsional Vibration. — ^The formula for a 
simple torsional vibration is derived in a manner exactly 
similar to that for a simple linear vibration given in § 10 
and is entirely analogous to it. 

For linear vibratory motion as in § 10, 


„ W Wd^-y , 

rorce = — a = — = — ky 

g g dt- 


(25) 


For the corresponding angular motion, such as that for 
the disk supported on a circular rod, performing a tor- 
sional vibration as shown in Fig. 14, 


Torque = L = -a = ~~ = - UQ . . (26) 

as in § 10. 

In the first case, a = linear acceleration and k = linear 
elastic constant. 

In the second case, L — torque, a = angular accelera- 
tion, and L-i = angular elastic constant = torque per 
radian of twist. 

Thus, for the angular oscillation, in correspondence to 
the linear of § 10, 

Frequency = / = .... (27) 

Ztt ' 1 


The quantity I corresponds to W of the first case and 
is recognized as the moment of inertia of the oscillating 

24 
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mass. In Fig. 14 it is the moment of inertia of the disk 
about its central axis. The moment of inertia of a body 
about a given axis may be thought of as its angular mass 
about that axis. 

24. Torsional Elastic Constant of Rod. — ^The torsional 
elastic constant of the circular rod of Fig. 14 is found as 
follows: 

Since Li is defined as the torque per radian of twist for 
the whole rod of length /, it is evident that the section of 



Fig. 14. — ^Torsional Vibration. Fig. 15. — Section of Twsted Rod. 


unit height shoum in Fig. 15 is twisted through 1// radians 
when L-i is applied. 

The elementary area rdddr at A in the side view of the 
unit section (Fig. 15) is sheared across the corresponding 
area at B by an amount determined by the product of the 
angle \/l and the distance r of this area from the neutral 
point at the cylinder axis. Since this is by definition the 
shearing strain e, (see § 15), 

r 


. ( 28 ) 
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Since by § 15, the shearing stress r = Ne„ where N is the 
modulus of shear, 


Nr 


(29) 


where r = the intensit>' of shear or the shear per unit area 
at the point where the elementary' area was taken. There- 
fore the actual shearing force on this elementary' area is 


rrdOdr = 


Kr-dOdr 

1 


• (30) 


The torque dLi produced by' this force is equal to the force 
times its distance r from the axis, or 


Therefore, 


dLx = 


Nr-dOdr 

I 


(31) 


2 , = /’*" f‘ Nt^dedr 
Jq Jq I 


-Na^ 

21 


. . (32) 


where a = radius of rod. 

This is the desired value of Li, the elastic constant of the 
rod. 

Substituting (32) in (27) it is found that the frequency 
of the system of Fig. 14 equals 


/ = 


1 J -AVg 

2 - ' 211 


C- hrNg 
2 -^ 211 


(33) 


25. Moment of Inertia. — The phy'sical meaning of the 
moment of inertia or angular mass of a rigid body may' be 
well shown as follows. 

Consider a rigid, weightless rod pivoted at 0, Fig. 16, 
wth a small weight 'jj on its end. A torque L about 0 pro- 
duces a force/ on a' such that 

^ w d^y 
■' g dP 

where y is measured along the path of w. 
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But 


^ . wr d^y 


. . (34) 


Also if 0 = angle swept out by r when w moves a dis- 
tance y, y = ?'0j so that 


d^y _ d^O 
dip ^ 


(35) 


Substituting (35) in (34) 

wr"^ d^d 


(36) 


Comparing (36) with (26) it is seen that the moment of 
inertia I of the simple system of Fig. 16 about 0 = wrP 
It is evident that an ^ ^ 

extended rigid body may Y 

be regarded as made up of <? r 

a large number of small 16 , — Simple Rotating System, 

weights w, each at its own 

distance r from the axis of rotation so that for the whole 
body I = sum of all these quantities. 

Expressed mathematically, 

j = = TFr" (37) 

W is the total weight of the whole body and r a particular 
value of r such that, if W were concentrated at a point, r 
would be the distance from W to the axis of rotation. The 
distance r is called the radius of gyration, and is seen 
to be the root mean square value of r determined by equa- 
tion (37), when the body is of constant density. 

In the exact process of evaluating I the summation indi- 
cated in equation (37) is performed by an integration, the 
component weights w being taken as infinitesimally small. 
The expression for I in this case is written 

^ =f f f = Wr^ 


(38) 
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where p is the density of the material and dV an element 
of volume. 

26. Routh’s Rule. — It is seen that when the radius of 
g^'ration r of the body about its a-xis of rotation is once 
known, I is easily found by multiplying r~ by TI’, the weight 
of the body. 

A verj' useful rule for finding r~ about an axis of sj'm- 
metiv' of an isotropic body, wthout an integration, is 
given bj’ Routh, being known as Routh's rule. It is stated 
as follows: 


sum of the squares of the perpendicular semi a.\es 
3, 4, or 5 


the denominator to be 3, 4, or 5 according as the body is a 
rectangle, ellipse (including circle), or ellipsoid (including 
sphere). 

For instance, for a rectangle with axis through 0 per- 
pendicular to the paper. Fig. 17, 



Fig. 17. — Reel' 
angle 


-2 0= + ft- 

r. = — r — 


(39) 


If the axis coincides with b, 
a- 

] ~ 3 


r = 


(40) 


where the thickness of the rectangle squared can be neg- 
lected in comparison with a~. 

For a cylinder of radius r about its axis. 


r = 




r- 

T 


. . (41) 


27. One Degree of Freedom with Damping. — Con- 
sider a system of one degree of freedom like that shown in 
Fig. IS, where, in addition to a spring of elastic constant k, 
a dash pot is introduced in parallel with the spring. The 
friction thus produced is assumed to be directly propor- 
tional to the velocity of motion. This law of friction is 
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taken because, besides being characteristic of dash pot 
action, this law permits of simple analytic treatment. 

If b signifies force per unit velocity, the force produced 
d'V d'V 

bv the velocity -r is The force produced by the 
^ dt at 

spring is ky. 

These two external forces which act on 
the weight are called the impressed forces. 

It is the impressed forces which cause the 
weight to accelerate in the direction in 
which they act. Both of the impressed 
forces are negative in this case since when 
the velocity and displacement are taken Fig. 18.— One De- 
positive they both act toward the origin Freedom 

. . . , . . ° with Damping. 

or point oi zero displacement, giving an 
acceleration towards the origin, 
ally, 



Expressed mathematic- 


d'^y dy . 


+ + = o . . . . (42) 


The symbol m = mass is used in place of W/g to sim- 
plify the notation. 

A solution of this equation is 

_iL 

y — As 2m sin (wo^ + 0i) . . . . (43) 

^2 

where wo = — and where A and 6i are arbitrary 

constants determined, respectively by the initial amplitude 
and phase angle when ^ = 0, 

Note that the vibration is harmonic and gradually 
dies out with time. If & = 0, expression (43) reduces to 
the form of equation (4), the vibration continuing indefi- 
nitely, and wo reduces to Wc = yfkjm. 
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28. Sustained Vibration with One Degree of Freedom 
and Damping. — For the case of a sustained vibration with 
one degree of freedom with damping, consider the case of 
Fig. IS with a periodic force Fcos a/ applied to the mass 
in an up-and-down direction. The stimulating force is an 
additional impressed force and has a positive sign, because 
it gives 7/1 a positive acceleration component. With this 
additional force, equation (42) becomes 

7/1 ^ -f -I- hy = F cos at . . . (44) 


The solution of this equation is 


y^Ae sin (W + <?,) + ^ _ (45) 


where 


tan /3 = 


ba 


k — 7/1 or 


V{k - war)2 4. 
VA/zik - *2 


and 0)0 = 


27/1 


and where A and 0\ are constants depending on initial 
amplitude and phase. 

The vibration is seen to be made up of two components, 
one on top of the other. The first is the free vibration 
component, whose amplitude and phase depend on how 
the vibration was started, and which rapidly dies out with 
time. The second is the forced or sustained vibration 
component which always has the same amplitude regard- 
less of the Initial conditions. 

As an example, consider the special case where 


dy 

y = 0 and -f- = 0 when / = 0 and w = 
at 



From these conditions, 
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Substituting these quantities in (45), 



F _M. . . , F sin Wct 

— — — 6 2m Sin coot H ; 

OCOo OcOc 


. . (46) 


This represents a vibration which gradually builds up to 
the value 


F sin ojci 



(47) 


which is reached when the transient term drops out. 
Since b is small in ordinary cases, coe and coo are very 
nearly equal to each other so that by the time the phase 
of the two components becomes opposite the amplitude 
of the transient has become almost zero. 



CHAPTER V 


SOLUTION OF VIBRATION PROBLEMS BY THE USE 
OF COMPLEX QUANTITIES WITH APPLICATIONS; 
VIBRATION OF COMPOUND SYSTEM 

2P. Use of Complex Notation. — The use of complex 
notation in solving for the steady state is well illustrated 
by the example of the previous section. Instead of repre- 
senting the exciting force b}' F cos oit it is taken as 
where j — V — 1 • Since imaginary quantities are of an 

entirely different nature from real, 

0 the usual convention in graphical 
jFsinat representation is to plot them 
g along the vertical axis and the real 
Fcos cut quantities along the horizontal axis 
so they are independent of each 
other. In Fig. 19, the , quantity 
Fig. 19.— Rc.iI and Imaginary F COS wl is given by the projection 
Components. of the revolving vector F on the real 

or horizontal axis, whereas the 
quantity Fe^"'' — F{cos wt -f j sin uf) includes the imaginaly' 
component as well as the real component of F, and thus 
represents the complete revolving vector. It is easily 
shown that 

= cos co/ -f- y sin co/ 

and ... (48) 

c " = cos (ot — j sin cot 

by expanding the separate terms into their corresponding 
scries, and comparing the resulting right- and left-hand 
members. If ut is measured in a counter-clockwise direc- 
tion, represents a unit vector revolving in the counter- 
clockwise direction and a unit vector revolving in 
the clockwise direction. 
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Equation (44) may thus be written as 

+ + . . . ( 49 ) 


regarding each term as a vector. Neglecting the transient 
state, the steady state of vibration may be assumed to be 
y = Ae^“', from which 

dy . , d^y 

and 


which shows the phase relations of the individual terms 
on the basis of the solution assumed. 

Substituting these quantities in (49), there is obtained 


(— jbo} + ^) y = 




. (50) 


Thus 


p . 

y = — where Z' = — mtsP- + k jboi 

Zj 

F 

The solution assumed is correct provided A = where 

Zj 


Z' is a complex quantity, so that A must be one also. 

The complex quantity Z' will be called the mechanical 
impedance of the vibrating system after A. G. Webster 
(reference 3). The mechanical impedance is that quan- 
tity by which F must be divided to give the vibration 
amplitude. The complex character of the amplitude factor 
merely means that it differs in phase from the exciting 
force F. 

Let us find the phase angle. This complex quantity A 
is of the form 


F 

a A- jb 


(51) 


where a is the real part and b the imaginary part of the 
mechanical impedance. Multiplying the numerator and 
denominator each by u — jb, there is obtained 


a - jb 


P ( 


a 


b 
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F 


Va- + b- 


(cos P — j sin /3) = 


Fc-^^ 

Va- + b- 


• ( 52 ) 


b , bct} 

where tan /3 = - is ■; for the problem just consid- 

a k — nio}- 


ered. Thus (51) represents a vector of magnitude 


— y which lags behind the vector F by an angle B 

V^F+lr- 

(see Fig. 19). Note that the solution of (44) given by 
equation (45) shows e.\actly this result for the steady 
state term as obtained without the use of comple.v quan- 
tities. 

30. Remarks on the Complex Quantity Method. — 
The comple.x quantity method, in mathematical language, 
is a shorthand process for solving for the particular integral 
of the linear diflferential equations of vibratorj’' motion 
where restoring forces are directly proportional to displace- 
ments. The motions of any vibrating system with lumped 
masses, stiffness factors, and friction coefficients can be 
completely expressed by such a set of simultaneous differ- 
ential equations, the preceding section containing an 
example for one degree of freedom. The same may 
be said of alternating-current electrical circuits with 
lumped inductances, capacities, and resistances, analo- 
gous respectively to masses, elastic coefficients (reciprocal), 
and friction constants of mechanical vibrating systems. The 
use of complex quantities in alternating-current electrical 
circuits is thus completely analogous to the process pre- 
sented in the preceding section. The alternating electric 
current amplitudes correspond to alternating velocity 
amplitudes of the masses of the mechanical system, and 
alternating charge displacement amplitudes correspond to 
mass displacement amplitudes of the mechanical system. 
The mechanical impedance is seen to differ from electrical 
impedance by the factor jw, because in mechanical systems 
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the impedance is so defined that force times impedance 
gives amplitude, whereas in electrical systems force times 
impedance gives current, which is the analogue of velocity. 
Referring to § 29 it is seen that sinusoidal amplitude and 
velocity differ by the factor yw, the latter being the deriva- 
tive of the amplitude, where j takes account of the 90° 
phase displacement between the two. 

Note that this process applies only where the steady 
state of vibration is required. When a mechanical system 
is acted upon by an alternating force of steady amplitude, 
in ordinary cases transient vibrations drop out in a few 
moments and are of no particular importance. The same 
is true for alternating-current circuits. Thus a process 
which neglects transients is sufficient. 

In solving mechanical systems it is best to proceed 
directly without the use of equivalent electric circuits, 
which are a superfluity, though often helpful to one versed 
in electrical theory. 

The impedance must be a comple.\ quantity in general 
because the periodic force and the resulting motions are 
vector quantities having different phase relations. The 
equations which are solved to obtain amplitudes are there- 
fore vector equations, the impedance being a comple.x 
multiplier which makes the two sides of each equation 
equal in both phase and magnitude. 

31. Rules for Solving Problems. — Draw a sketch. of the 
mechanical system, as shown, for e.xample, in Fig. 22, where 
the dash-pot indicates damping. Specify the amplitude 
of motion of each junction of elements by yi, y 2 , ys, etc. 
A mass has only one point of application while springs and 
dash-pots have two. The intersections of all springs and 
dash-pots attached to a mass are identical. A spring 
and a dash-pot may be connected without any mass 
at the junction. Note from § 29 that the mechanical 
impedance of a single mass m is — moy^, that of a single 
elastic connection = k, and that of a single damping factor 
(friction proportional to velocity) = jbo}. Set up a set of 
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force equations, one for each junction of elements according 
to the following rules: 

Rule 1. Multiply the amplitude 3 ' of the junction by 

everj’’ impedance element which is connected to it. 

Rule 2. Subtract from this quantity the product of 
the amplitude of the farther terminal of each 
spring or dash-pot into their respective impe- 
dances. 

Rule 3. Set this quantity equal to the force applied 
at the intersection if it is present. If more than 
one force Is applied to the system, account must be 
taken of their phase relation. 

A set of linear equations in 3 * 1 , yo, etc., will be obtained 
whose number equals the number of individual junctions 
of the system, which can be readily solved for amplitudes 
by determinants. 

Rule 4. To find the phase angle of the displace- 
ment with respect to the exciting force, if 
required, put the expression for that displacement 
in the form y = F/{A + jB) which always can 
be done. Then, as outlined in § 29, y has the 
magnitude Ffy/ A- -f and lags behind F by 
the angle /3 where tan /3 = B/A. 

Any problem of this type can be solved by this process, 
and most practical problems with one degree of freedom 
will be found to solve very easily (reference 4). 

Note that this analysis applies to motions which take 
place onlj'^ in one independent coordinate (y in this case). 
Motion in other coordinates is treated separately when 
required. 

32. Illustrations of Set-up of Equations. 

1. Periodic force F cos (at acting on a mass m, Fig. 20. 

Let amplitude of m — yj. Then 

— F 

- ymor = F, yi = — - . 

moP- 


. . ( 53 ) 
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This equation gives the amplitude yi with a minus sign 
which means that force and displacement are in phase oppo- 
sition. If a weight held in the hand is rapidly oscillated to 
the right and left this phase opposition is at once perceived 
in that the maximum right-hand force is exerted when the 
weight is in the extreme left position, and vice versa. 

2. Periodic force F cos oit acts on spring, Fig. 21. 

In this case 

F 

yxk = F] yy 


Here the force and displacement are in phase, as a little 
consideration will show must be true. 


I Fcos wt 

displacement 


I Fcos lot 

= Maximum 
displacement 


c=>< 



Fig. 20.— Vibration of Single Mass. Fig. 21. — Vibration of Single Spring. 


3. Take the case of Fig. 18, where a periodic force F 
acts upon a mass m supported by a spring k, with a dash pot 
of friction coefficient b as shown. Since there is one moving 
junction point, one equation is sufficient. By rule 1, we 
have 

yi{k — mop -\-jbui) = F 


No further terms are required because the lower ends of k 
and b are fixed. (See rule 2.) 

By rule 3 the left-hand member is set equal to F. 
Solving for yi 


yi 


By rule 4 


F 

k — moP jbo) 


yi 


V (fe — moP)^ -}- b'^oP 


where tan /3 


. . . (54) 

b 

k — moP 
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and where means the phase angle is measured back- 
wards against time, that is, Vi lags with respect to F. This 
problem is seen to be the same as the one treated in § 29. 

4. Take the case of the previous problem e.\cept that 
k and b are in series as shown in Fig. 22. 

In this case there are two moving junction points, one 
at in whose amplitude of motion is specified by ji and one 
between the spring and dash-pot whose motion is specified 

by >’ 2 ._ 

Using rules 1, 2 and 3, the following pair of simultaneous 
equations is obtained : 





w 

T 


Vi 


y, 


yi(jbo} — 7;; 03-) — y^jbo} = F 

y^ik jbu) — yijbu = 0. 

From the solution of these 
equations the displacement 
amplitudes yi and y 2 can be 
found, and by applying rule 4 
their phase angle with respect 
to F is found. 

33. Compoimd System. — 

Fig. 22.— Spring Fig. 23. —Com- 'j'jjg amplitudes of vibration 
and Dash-pot in pound System. r . i .1 . . r 

for the steady state of tlie com- 
pound system of two degrees 
of freedom shown in Fig. 23 are easily found by the use of 
comple.\ quantities. 

The general procedure in a case of this sort is as 
follows. 

Use a separate symbol for the amplitude of motion of 
each mass. In this case ji is the amplitude of nii and y 2 
that of nio. 

Set up a separate equation for the forces on each mass. 

In this case, bj'^ rules 1, 2, and 3, § 31 


yi(— miur -f ki) — 3'2^i = F| 

3-( — in2oP -f- “F ^ 2 ) — yiki = Oj 


( 55 ) 
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In this problem no damping is present, so that the 
symbol i does not enter. 

By rule 1, the amplitude of each mass is multiplied by 
each impedance which is connected with it. For the mass mi 
this gives yi(—Wiw“ +^i). Then from this quantity is sub- 
tracted the product of each of these impedance compo- 
nents into the amplitude of their other terminals. For 
mi there is only one such term, — y 2 ki. Set the sum of 
these terms equal to the exciting force F which acts on nii. 
This accounts for the first equation. The second is obtained 
by the same process. Three impedances are seen to act at 
W 2 . Only the spring of elastic constant ki has its farther ter- 
minal move a finite distance, which terminal moves through 
the amplitude yi so that the term yiki is subtracted in this 
case. 

These rules may easily be established by setting up 
simultaneous differential equations in the usual way and 
making substitutions as in equation (50) and rearranging 
terms. 

The amplitudes are most easily found by solving the 
simultaneous equations by the use of determinants. For 
this case, first arrange the equations in symmetrical form, 
using ai, bi, G 2 , and 62 as coefficients. 


(iiyi + a 23'2 = F 
biyi -f- h2y2 — 0. 


Then solve for yi and y 2 in the usual way. 


3'i = 


F 

(I2 


0 

b2 

Fb2 

Cl 

bi 

a2 

b2 

C162 — C261 

Cl 

F 


bi 

0 

^ -Fbi 

Cl 

bi 

02 

b2 

aib2 — C2J1 


. (56) 


( 57 ) 


3-2 = 
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Substituting, 


F { — + k\ + A’2) 

(— vi\(jr + ki){-— J«2W- + ki +A’2) — A’r 

^ 

(— wiur + ^»2W- + ki + ko) — kr 


( 58 ) 


The resonant frequencies of this system for which both 
3 'i and 3'2 become infinite are found by setting the denom- 
inator equal to zero. This gives 


/= 


l-y/l mi mi ' W2 mi mi) mmi 


(59) 



Fig. 24. — Compound S>-s- 
tcm with Damping. 


The system thus has two resonant 
frequencies, for the lower of which 
the two weights vibrate in phase and 
for the higher of which they vibrate 
in phase opposition. This is seen to 
be the same system as that discussed 
in § 12. 

34. Compoimd System with Damp- 
ing. — This system is shovTi in Fig. 24. 
Using the general rules of § 31, the 
force equations corresponding to (55) 
become 


yi(-?nic<r+/:i-h7&iw) —yi{ki+jbi(S) = F 

yi(^ — 7n2W‘ "hAi “}“A2~1 'J&i w-Fj& 2 w) — yi {ki ~\~JbicS) = 0 


(60) 


The solution of these equations is carried out in e.\actly 
the same way as outlined in the previous section. 

With a little practice the equations can be readily set 
up for any system of this t 5 ’-pe, however comple.x. The 
actual e.xpressions for amplitude may be very complicated 
indeed. Nevertheless this is the simplest and most direct 
method knoum for handling such problems. 



CHAPTER VI 

LONGITUDINAL AND FLEXURAL VIBRATIONS OF RODS 


35. Longitudinal Vibration of Rods. — Before taking up 
flexural vibrations of rods which are important in vibration 
work the interesting case of longitudinal vibrations of rods 
will be treated, which is somewhat simpler than the former. 

If a tension is applied in the X direction on the piece 
of rod shown in Fig. 25 fixed at the origin the section AB 
will be displaced a small distance to the position CD, 
owing to the elastic strain. This displacement is called 
the total strain at AB 
and will be designated 
by the symbol u. It is 
important to appreciate 
that this differs from 
the strain per unit length 
e which determines the stress intensity 
definition that 




i 




.BD 

Fig. 25. — Longitudinal Strain in Rod. 


It is evident from 


otram = € = — 


dn 

dx 


(61) 


dx 


SEe 


SE(e*d€) 

>- 

de 


Fig. 26. — Section of Rod. 


to the 
strain e 


Let us find what forces act 
on the element or section of the 
rod of thiclcness dx shown in 
Fig. 26 enlarged. If ^ = cross- 
sectional area of the section and 
the strain e increases in the 
positive direction of X the force 
on each side of the section due 
is as shown in the figure. The 
positive so that the stress is a 
41 


elastic strain 
is assumed 
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tension. The resultant axternal force acting on tliis 
section is 

SE(^i + 

If p = mass per unit volume of the rod material, the mass 
of the slice of Fig. 26 = pSdx. Therefore 

oj err 9~R , 

pSdx — V = SE — dx = SE — - dx 

^ 3/- dx 3.V- 


This is the mathematical statement that the mass of the 
slice times its acceleration is equal to the resultant external 
force on it. 

Rearranging and cancelling common factors, 


P 


d-n 

di- 




. (63) 


Although here derived for a bar fixed at one end, the 
general solution of this differential equation gives all pos- 
sible kinds of longitudinal elastic waves in such a rod, 
regardless of end conditions. 

This vibration system differs from any thus far con- 
sidered in that, instead of separate masses joined by elastic 
connections, a continuous elastic medium of uniformly 
distributed mass is involved. The vibration character- 
istics of such a medium are the same as though it con- 
sisted of a verj' large number of very small separate masses 
all connected by elastic bonds of equal strength. The 
continuous medium is the limiting case where the number 
of the component masses is indefinitely great and their 
size indefiniteh' small. A system with an infinite number 
of elastically connected masses should have an infinite 
number of degrees of freedom. It is very interesting that 
this is just what the solution of equation (63) shows. 

36. Solution of Equation for Longitudinal Vibrations. — 
The form of solution of equation (63) which satisfies the 
requirements of this problem is 

u = . . . (64) 
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where Ai and A 2 and depend upon the end conditions of 
the rod, and w, the frequenc}’- constant = 27r/, depends 
upon its physical properties. This solution is given di- 
rectly witliout going into the method of obtaining it, which 
can be found in books devoted to differential equations. 

For a free-free bar (both ends free), since the stress 
Ee must be zero at the ends, 

dtt 

e = = 0 when re = 0 and when re = / . (65) 

dx 

where I = length of the rod. 

— = = 0 when re = 0 

dx 

Therefore, Ai = A 2 = A /I, say. 

Then from (64), 

u = d -I- .... (66) 

From equations (48) 

u = i4e^“‘ cos iSre = i4(cos cos jSre -f j sin co/ cos |3.x'). (67) 

Equation (67) gives the amplitude of vibration n as a 
vector revolving at a constant angular velocity co whose 
length varies along the rod according to the law cos jSre. 
(See §29.) 

Since the actual vibration amplitude is the projection 
of this vector on the axis of real quantities, the 7 compo- 
nent of (67) can be dropped. Its use is of value merely 
in simplifying the process of solution of the differential 
equation. In specifying the amplitude u it has no physical 
meaning. When thej component is dropped, (67) becomes 

u = A cos w/ cos ^x (68) 

When (68) is plotted it is found to represent a series 
of standing waves along the rod, as showm in Fig. 27, 
where the amplitude u is plotted along the Y axis. These 
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waves rise and fall between the maximum and minimum 
positions shown in the figure with fixed nodal points. 

If no = amplitude when / = 0 at .r = 0, from (6S) A = Uo, 
so that 

11 = Uo cos u! cos /3.V .... (69) 

Also 

du „ . • „ 

— = — Uop cos ut sm px 
ax 

which, from (65) when .v = I, 

= — Uo^ cos cot sin jS/ = 0 

This requires that sin fit = 0, since cos cot is not always 
zero. Therefore 

= n-ir; ^ (70) 

where n has any integer value. 

Thus from (69) when .v = /, u = Uo, as indicated in 
Fig. 27. 

Since /3 may have any one of an infinite number of 
values the rod may vibrate in any one of an infinite number 

of modes of which Fig. 
27 represents the partic- 
X ular case where n = 4. 
Of course, the general 

Fro. 27.— Vibration in Four Nodes. solution (69) is the SUm 

of all harmonics cor- 
responding to all values of n. These may, as brought out 
in § 12, be superimposed upon each other in any combina- 
tion without their individual frequencies being alTected. 

From Fig. 27 a complete wave length, designated by 
X, is equal to twice the distance between nodes. It equals 
the value of .v corresponding to one cycle of cos ^x, that is 

when /3.V = 2;r; x = X = ^ = -^ . , (71) 
since from (70), /3 = 

I' 
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Thus for Fig. 27 where n = ^,\ = 112. 

The only constant left to evaluate is the frequency 
j _ ^/27r. This is evaluated by substituting the solution 
(69) in the original differential equation (63), which it 
must satisfy. This gives the relation 

— pno(jy^ cos ul cos cos ut cos = 0 


from which 


K mr [E 



CO = ^V- 
/ 

7 T V 


Therefore 

/■I'/ 

E 

• • • • • 

p 

. . . (72) 


It is seen from (72) that each harmonic has a different 
frequency, proportional to the number of the harmonic or 
to the number of fixed nodes. For the case of Fig. 27 
where n = 4, 




Now, for a progressive train of waves, evidently fre- 
quency equals velocity of motion of the wave train divided 
by wave length, or expressed mathematically 



where V = velocity of the "wave train. 
From (71) and (72) 



This equation gives the velocity of transmission of a wave 
train of longitudinal vibrations along the rod. It is seen 
to depend only upon the elastic modulus E and the density 
P of the material of the rod. It is really the velocity of 
transmission of sound along the rod, which for steel is 
equal to about 17,000 feet per second, or about 16 times 
that of sound in air. 
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37, Travelling and Standing Waves. — It will be here 
shown that a train of standing waves such as that obtained 
in the solution of the problem of the previous section results 
from the superimposition of two trains of travelling waves 
of equal wave length moving in opposite directions (reflect- 
ing back and forth along the rod). 

Let Fig. 2S represent part of a moving wave train here 
shown when one of the travelling nodes is passing through 
the origin. 

If the wave were fixed in this position its equation 
would evidently be 



2 

Fig. 28. — ^Travelling Wave. 


It is seen that when the distance .v is a multiple of the 
wave length X, the wave angle is a multiple of It:, as it 
should be. 

In order that this wave train shall move to the right in 
the positive .v direction the wave angle corresponding to 
any point on the .v axis such as the origin must continu- 
ously decrease so that more and more negative portions of 
the war-e progressively pass it. If this wave velocity is 
V, this decrease of wave angle per second must be Vt/'K, 
so that the equation for a Avave moving to the right Avith 
the A'elocity V becomes 

y =^siny (.V - FO .... (75) 

In like manner, if the AA^aA^e train is moving to the left, the 
AvaA'e angle corresponding to each point .v must progres- 
siA'ely increase so that the equation for this AA'^ave train is 

y = Y sin ^(x -f FI) .... (76) 
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To find the type of wave produced by superimposing these 
two wave trains, add (75) and (76) after putting in the 
form for the difference and sum of two angles. The 
resulting equation is 

. Ittx 2TrVl 

y = yosm^— cos — — . . . . (//) 

A A 

This equation is seen to be identical in form to (69), 
which has already been shown to be the equation of a 
standing wave train. Since, from the previous section, 
\ = = V/f, making these substitutions, (77) becomes 

y = yo sin cos ut (78) 

It is often very useful in vibration work to separate such 
a standing vibrating wave train into its two oppositely 
moving wave trains. The analysis just given furnishes the 
basis for this. 

Note also that a continuous and complete reflection 
of these wave trains from the ends of the bar is necessary 
for this agreement. 

In the previous section the velocity of wave motion 
V = s/EIp was obtained from a standing wmve train. 
When it is appreciated that this velocity is the actual ve- 
locity of the two oppositely moving wave trains of which the 
standing wave train is composed, this velocity assumes 
greater physical reality. 

The expression for the velocity of sound in air is the 
same as (73), where E is the elasticity of the air. 

The isothermal elasticity of air equals the pressure p. 

Since in sound waves the compressions and rarefac- 
tions are relatively instantaneous, the adiabatic elasticity 
must be used = jp where 7 = Cp/c^, the ratio of specific 
heat of constant pressure to specific heat of constant 
volume. For air, formula (73) becomes 

V = 

^ p 


(79) 
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If mass per unit volume w/g is substituted in place of 
p, equation (79) becomes 

= (80) 

Another point of interest is that when a wave of com- 
pression travels along a rod of rigid material like steel it 
is attended by a wave of lateral expansion of small mag- 
nitude because of Poisson’s ratio. This means that waves 
of sound in an extensive rigid medium travel at a slightly 
higher velocity than that given by (73) because in such 
a case no lateral expansion is possible, thus giving a 
slightly higher effective value of E. The analysis of this 
case can be found in works on elasticity (reference 5). 

38. Flexural Vibrations of Rods. — Let Fig. 29 represent 
part of a rod or bar subject to transverse forces in the T 


Y fdx 



t F QM 

1 


1 


^0 

t 


1 

(l 




Fig. 29. — Forces in FJc.xural Vibration. 

direction, the axis of which coincides with the X axis when 
it is undeflected. 

Let / = transverse force per unit length, positive 
upwards; F the shearing force on any section, positive as 
shown; and M the bending moment on each section,’ 
positive when acting clockwise on a section as viewed in 
the positive X direction. 

The shear on any section evidently equals the integral of 
Jdx from the origin to that section plus the shear Fo at the 
origin as seen from the figure. Expressed mathematically, 

F = ^Sdx + Fo .... (81) 
Also the change in bending moment M per unit length at 
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any point equals the shear at that point. This can be 
seen to be true because unit length has unit moment arm. 
Expressed mathematically, 


dM 

dx 


= F 


(82) 


From (81), dF/dx = f, so that 

d^M 

dx^ 

Now since, from (23), 


(83) 


M = El 


d^y 

dx^ 


dx^ dx^ ^ 


. . (84) 


Thus the fourth derivative of the deflection curve at any 
point is equal to the load per unit length at that point. 
The relation (84) will be made use of in analyzing fle.\ural 
vibrations of rods. 

Assume the rod to be performing sinusoidal or har- 
monic flexural vibrations. 

Thus the motion of the rod at some point x is equal 
to y = u cos wt, where u = the maximum value of y or 
the amplitude of the motion, li w — weight per unit 
length of the rod the weight of an element of length dx 
equals wdx. Its acceleration due to the sinusoidal motion 
is 


IF 


= — U(j}^ cos (at 


The force producing this acceleration of the weight wdx is 


wdx d^y 
g dfi 


wdx „ 

Uca^ cos (at 


( 85 ) 


g 
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If /' = the accelerating force per unit length at that 
point, when y = ic 

7iu}~ = / ax 

g 

or 


w 


- - = /' . . . 
g 


. (S6) 


Notice from Fig. 30 that/' acts downwards in accelerating 
the element w towards the X a.\is. The force / is caused 


fdx 
f /uidx 

mit:— 


-t- 

fdx 


Fig. 30. — Action and Reaction of Element. 


by the reaction of w against the rod structure and is thus 
equal and opposite to /'. Thus 


j«co2=/ (87) 


The differential equation from which the frequency con- 
stant w can be found is obtained by substituting (87) in 
(84) and setting y equal to «, the maximum value of y for 
which (87) was derived. This gives 


Or 


where 


d'a 

dx^ 


wco- 
-777 a 
gEI 


0 


( 88 ) 


d'n 

— viHi = 0 
dx- 


(89) 



(90) 


39. Solution of Equation for Flexural Vibrations. — The 
general solution of (89) is well known to be 

u = A cosh 7}ix -j- B sinh mx C cos mx + D sin 77ix. (91) 
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Nole that 


cosh .%* = ^ ^ — and sinh .t = ^ n — (^2) 

For la bar fixed at one end and free at the other, called a 
fixed-free bar, 

when .r = 0, = 0 and — = U 


d-u dhi 

wlien:c = i. ^ = 0 and ^ = 0] 


. (93) 


since at the free end of the rod the bending moment 
M - Q and the shear = dM/dx must also be zero. 

Applying these four end conditions to (91), four sepa- 
rate equations of condition, as they are called, are obtained 
from which the four arbitrary constants A, B, C, and D 
can be eliminated, resulting in the final equation. 

cosh ml cos ml \ = 0 . . . . (94) 

This equation is satisfied by an infinite number of values 
of ml, each of which corresponds to a particular mode of 
vibration. Since I is 
known, m can be 
found. The fre- 
quency of that mode 
can then be calcu- 
lated from (90) . The 
shape of the deflection 
curve at maximum Fig. 31. — Modes of Vibration of Cantilever, 
amplitude can be 

found from (91) when the constants A, B, C, and D have 
been evaluated. (See § 88.) 

Figure 31 shows the first three modes of vibration with 
the corresponding values of the constant ml. 

First mode, ml = 1.875 
Second mode, ml = 4.695 
Third mode, ml = 7.856 
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For higher modes, vil = (2« — l)(7r/2), to a close app ;oxi- 
mation, n being the number of the harmonic. 

For a supported-supported bar, the end conditions are 

« = 0 and = M = 0 when .v = 0 or x = I (95) 


These conditions applied as before yield the simple relation 

sin ml = 0 (96) 

This gives ml = vx, where n = number of the mode. 



n = I 


7J. = 2 

71 = 3 


ml = ff 


ml = 2ff 


ml = 3 tX 


Fig. 32. — Modes of Vibration of Supported Rod. 


Figure 32 shows the first three modes. 

n = 1 ml = TT 

« = 2 vil = 27r 

n = 3 ml = 3 t 

It will be shown later that this analysis holds for the 
critical speeds of a uniform shaft. 

Note that in this analysis the cross-section of the bar 
is taken as uniform so that the second moment of cross- 
section 7 is a constant. If the rod is not uniform I will 
vary with .v and the solution may be very complex or 
impossible. 

For good treatments of the subject of vibrating bars 
see reference 6. 
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THEOREMS ON THE ROTATION AND TRANSLATION 
OF RIGID BODIES AND APPLICATIONS 

40. Theorems in Rigid Dynamics — Torque . — A pure 
torque or couple acting on a rigid body can produce rotation 
about its center of gravity only, regardless of where the couple 
acts on the body. 

Let a couple act on the body shown in Fig. 33 about an 
axis parallel to the vector A. Assume that the Y axis 
coincides with the axis of 
rotation of the body. As 
the body is accelerated 
about the F axis the ele- 
ment dxdydz is acted upon by 
the force df = {w / ^adxdydz, 
where w = weight of unit 
volume and a = linear ac- 
celeration of the element in 
the direction Z. 

If a = angular accelera- 
tion, a = ax. This is true 
whether the element be in 
the XY plane or not, as can 
be seen from the top view 
shown in the figure. Then 

*WOl. 

df = — xdxdydz (97) 

Since for a pure torque no resultant translating force 
exists, the forward accelerating forces parallel to Z on 

53 
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one side of the Y axis must equal the backward accelerat- 
ing forces on the other side, so that 

— J" J' J' xdxdydz = — J' J' J' xdxdydz . (98) 

{X positive) (.Y nc^niUvc) 

But this is the condition that the center of gravity lie in 
the YZ plane, namely, that the sum of everj’’ element mul- 
tiplied by its perpendicular distance from one side of tlie 
plane equal the same sum for the other side. 

By exactly similar reasoning the center of gravity may 
be shown to lie also in the XY plane, so that it must there- 
fore lie on the Y axis. Since by hypothesis the Y axis 
is the axis of rotation, the center of gravity must lie on the 
axis of rotation, and vice versa. 

41, Single Force . — A single force acting on a rigid body 
can cattse pure translation Toithout rotation only when its 
line of action passes through the center of gravity of the body. 

Assume a system of coordinate axes to pass through 
the body as in Fig. 33, but in this case take the origin at 
the point through which the force must pass to cause 
pure translation and let the force act along the Z axis. 
Ever}' element of the body is then subject to a linear 
acceleration a. In this case the reaction of each element 
causes a turning moment about the Y axis given by 

dM — "^xdxdydz (99) 

For pure translation the sum of the moments produced 
by elements on each side of the Y axis must be equal and 
opposite so that 

S J' ~ ~ J' J' J' xdxdydz . (100) 

(.Y positive) (A’ negative) 

As in the previous section the center of gravity must 
lie in the YZ plane. In like manner it is also found to lie 
in the XZ plane, which two conditions require that it lie 
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on the Z axis. If now the force is applied along the Y 
axis a third condition is obtained which proves the origin 
to be the center of gravity of the body. 

Since by hypothesis the origin was taken as the point 
through Avhich the force must act to cause no translation, 
this point and the center of gravity are coincident, and 
vice versa. 

Conversely a single force will cause a rigid body to 
rotate in all cases when its line of action does not pass 
through the center of gravity. 

It is evident that if a force must act in line with the 
center of gravity of a body to cause linear acceleration of 
the body vdthout rotation, the reaction to 
that force must also lie on a line through the /p 

center of gra\’ity of the body. 

Note from Fig. 34 that a single force F / [ 

acting on a rigid body at P, say, along a line ( 
not passing through its center of gravity is 
equivalent to a couple consisting of a pair of /p 'P 
forces F and moment arm AG together with 
a single translating force F passing through gi Force Acting 
its center of gravity. If the body is placed on Mass, 
in the figure so that the original force and G 
lie in the plane of the paper, the other two forces must 
also lie in that plane. 

It follows from this that, no matter where a single 
force acts on a rigid body, it always causes a translation of 
the center of gravity of the body as though the force were 
acting at the center of gravity. 

Rotation about the center of gravity is produced in 
addition because of a torque equal to the force times the 
distance of its line. of action from the center of gravity. 

A pure torque, however, can cause rotation only about 
the center of gravity Avith no additional translation, regard- 
less of where the torque acts on the body. (See § 40.) 

42. Further Remarks on the Motion of Rigid Bodies.— 
It is well known that a spinning body rotates about an 
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axis through its center of gravity. It is readily shown by 
a method similar to that used in § 40 that this is necessary' 
if the radial centrifugal forces are to be in equilibrium. 

It should be pointed out that, although the axis of 
rotation of a rigid body passes through its center of gravity 
when acted upon by a couple, the a.xis of rotation and the 
axis of the couple are not necessarily parallel. For instance, 
a couple represented by A in Fig. 35 acting on the end of a 
free bar will cause it to rotate about the axis represented 
by B. Every rigid body has three principal axes of inertia 
to one of which the couple must be parallel' in order that 
the axis of rotation and the axis of the couple shall be 
parallel. For the rod shown, one principal axis is parallel 



Fig.. 35. — Torque Axis and Fig. 36. — Equimomental Skelc- 
Rotational Axis _not Parallel. ton. 


to the axis of the rod. The other two are perpendicular 
to it. 

It can be shown that every rigid body can be repre- 
sented by an equivalent skeleton of three rods, as shown 
in Fig. 36, so far as its inertia characteristics are concerned, 
called an equimomental skeleton. The three principal 
axes of inertia coincide with these three rods, the greatest 
inertia in this case being that about C, called Jc, and the 
least Ia, with In intermediate. If all three of these happen 
to be equal, as for a sphere or cube, the axis of rotation 
and the axis of an applied couple always are parallel. 

For an excellent discussion of the principles presented 
in this chapter see reference 7. 

43. Illustrations from Common Observation. — com- 
plete discussion of these phenomena would take us far into 
rigid dynamics and can be found in works devoted to this 
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subject. For our purposes the outline in the preceding 
three sections is sufficient. The principles given are funda- 
mental and not difficult to remember when once clearly 
seen. They are illustrated by many familiar phenomena. 
We know that when a body is dropped it falls in a straight 
line with no rotation. By definition, gravity acts through 
the center of gravity of a body. Thus the truth of the law 
given at the beginning of § 41 is at once seen. 

A rotor in so-called dynamic unbalance will tend to 
wobble about its center of gravity owing to the centrifugal 
couple although its actual motion may be much modified 
by the elastic action of the bearing supports. This couple 
can be counteracted by a pair of balance weights producing 
an equal and opposite couple acting anywhere on the rotor, 
even out on the end of the shaft if there is room and if the 
shaft does not deflect appreciably because of the couple. 
This follows from the law given at the beginning of § 40. If 
a stick is thrown and rotated at the same time it will be seen 
to revolve about its center of gravity as nearly as can be 
observed, the center of gravity describing a smooth curve as 
though all the mass of the body were concentrated there. 

If a button is spun on a string at high speed it will 
revolve about its center of gravity even if it is somewhat 
eccentric on the string. 

If an unbalanced rotor on a flexible shaft is revolved at 
a speed far above its critical, its tendency to revolve about 
its center of gravity is so strong that the latter point re- 
mains practically stationary and the shaft is forced to 
revolve around it in spite of the elastic reaction of the shaft. 
The phase reversal which takes place in passing through 
critical speed is a result of the tendency for the center of 
gravity of the revolving system to remain stationary assert- 
ing itself against other forces. (See § 42.) 

In vibration work, occasion sometimes arises to find 
the magnitude of the periodic force which may cause a 
rigid body to vibrate, such as the pump mechanism of an 
electric refrigerator. Suspension by long, flexible springs 
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will permit the body to vibrate as though practically free 
in space, when its vibratory motion can be accurately 
obser\-ed. If the body is not symmetrical it may be found 
to perform torsional vibrations about an axis not parallel 
to the axis of a periodic exciting couple. This, of course, 
is because the axis of the exciting couple is not parallel 
to a principal axis of the body. 

Large errors have been made in estimating the magni- 
tude of the stimulus causing a vibration, such as the 
amount of unbalance in a motor, by not suspending prop- 
erly. If the motor Is placed on a pad on an ordinarj' 
table, the latter may amplify the vibration amplitude 
several times, and vice versa. This subject is discussed 
further under the heading of elastic suspension. (See 
§67.) 

44. Application to Problems. — It follows from §§ 40-41 
that the acceleration of the center of gravity of a rigid body 
is the same as if the translating forces acted at that point, 
and rotational acceleration takes place about the center of 
gra\'ity as if the body were pivoted at that point, regard- 
less of the actual translation which may be present. 
Thus, to find the motion of a rigid body acted upon by a 
gi^'en set of forces, translation and rotation may be treated 
separately. 

First: set up force equations for acceleration of the 
center of gravity. 

Second: set up equations for torque about that 
point. 

In general, if all three dimensions are involved, there 
will be three equations of translation and three of torque, 
each of which may be independently set up. 

Many practical problems can be set up so that there 
is one axis of translation and one of torque. If the axis 
of torque is parallel to a principal axis of inertia the rota- 
tion axis is parallel to the torque axis, and the problem 
is comparatively simple. Otherwise the torque must be 
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resolved into components parallel to the principal axes of 
inertia, so that there are two or three torque equations 
instead of one. Otherwise the procedure is the same. 
A torque can be represented by a vector parallel to its 
axis and resolved as though it were a single force. 


Example 1. — Find the acceleration 
plane inclined at an angle /? (Fig. 37). 

Leto = constant linear accelera- 
tion; 

a = constant angular accele- 
ration ; 

r = radius of cylinder; 

W = weight of cylinder; 

Jo = moment of inertia of cyl- 
inder about its axis. 


of a cylinder rolling down a 



Fig. 37. — Cylinder Rolling 
Down Incline. 


By Routh’s rule (§ 26) 


I = 


Wr^ 

2 


Let F = tangential force reaction on the cylinder from the 
plane. If a and a are positive for right-hand motion, 


for translation 


and for rotation 


W 

— a = W sin /8 


- F 


Wr2 


2g 


a — Fr 


Since a = a/r, from (102) 
Substituting in (101), 


F = — a 
2g 


( 101 ) 

( 102 ) 


31F , . 2 

a = PF sin |3. Therefore a = - ? sin 6 
2g 3 


(103) 


Example 2— A cylinder of radius n rolls back and forth on a 
concave surface of radius of curvature r 2 . Find the frequency 
for an oscillation of small amplitude. 
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Using the sjnnbols of previous problem with those of Fig. 38, 
Translation 



since s is measured from AB, and 0 is small, so sin 0 is set equal 
to 0. 

Rotation lFri“ dr4> 


= - Fn . 


. (105) 


2g dfi 

since 4> is measured from the middle position of the cylinder and 
the torque produced by F is thus negative. 

Let us express 4> and 6 in terms of 5. From the figure, 

5 


e = 

To — ri 

Also (<;> + 0)/r2 = 6/ri, from which 

r2 — r, s 

4> = 0 = — 

n ri 

Substituting (107) in (105) and solving for F, 
_ 

2gdfi • • • • 


Substituting (108) and (106) in (104), 

= 0 . 


3 d^s s 
2g dt~ r2 — ri 


(106) 


(107) 


(108) 


(109) 


This represents a simple linear vibration as given by equation (3) 
whose frequency is 


1 




2g 


. . ( 110 ) 


'3(r2 - r,) 

Example 3a. — Assume a cylinder to be removed from the 
earth so that it floats in space. How will it vibrate if a sinusoidal 
torque of ma.\imum value L be applied to it parallel to its axis? 

Translation = 0 


Rotation 0o 


{=?)■ 


Oo =- 


/co2 


. . (Ill) 
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where 0o = amplitude of the angular vibration about the axis of 
the cylinder. 

The torque equation is written down directly, using the 
mechanical impedance — {I/g)o3^ for rotational oscillations in 
place of — wco^ used for linear vibrations. 

(See § 31.) 

The sign of 6o is negative with re- 
spect to the torque because the two are 
in phase opposition, as can be easily 
shown by experiment. 

Example 3b . — ^Take the same prob- 
lem as the previous except that the 
cylinder is restrained by a pivot at P 
with a tangential spring support as 
shown in Fig. 39, whose total elastic 
constant equals k. Otherwise the cyl- 
inder is free. Find the natural or res- 
onant frequency. This problem differs 
from previous ones in that force of Cylin- 

translation and also the torque each 

depend upon both the linear displacement x and the angular 
displacement d, as shown in Fig. 40. 

The point P is displaced by an amount re -1- rd from its rest 
position Po. The tangential force (taken horizontal since 9 is 
small) is then equal to ^(x -f- r6). 




Fig. 39. — Oscillating Cylinder with 
Single Elastic Constraint. 



Fig. 40. Same as Fig. 39 Showing 
Displacement. 


Translation 

(x + r9)k — 

o 

II 

^ . . (112) 

Rotation 



(x-[-re\^ „ 

/Wr^ \ 


V ^ )^r^-9 

) = • 

. . (113) 
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Note that in rotation the elastic constant L\ = kr~, since it is 
defined as torque per radian. The angular elastic mechanical 
impedance is thus kr-, and the linear as before is simply k. Also 
the angular mass impedance is — {I/g)ur and linear mass im- 
pedance is (— Tr/g)£o-. 

In equation (113), since kr- is angular impedance, the dis- 
placement is expressed as an angle. 

As before, the use of impedance requires that these be vector 
equations. Since no j terms are present the vectors are either 
in phase or differ in phase by 180°. They must therefore be along 
the same line. 

Solving these equations by the usual method will give x and 0. 

Solving for 0, 



The resonant frequency is that for which 6 becomes infinite (and 
also x). This occurs when the denominator of (114) equals 
zero. Thus 


TTr2 „ 

C0“ 

2g 



from which 


^ W 


1 . 


(115) 


This system has one natural frequency and thus one degree of 
freedom, which is evident from inspection. 



CHAPTER VIII 


CRITICAL SPEEDS OF ROTORS— RAYLEIGH’S AND 
MAXWELL’S THEOREMS 


45. Critical Speeds of Rotors. — Consider a weightless 
shaft carrying a flysvheel of weight W and supported in 
bearings at A and B (Fig. 41). The up-and-down vibra- 
tion frequency is calculated as before for the case of a sys- 
tem of one degree of freedom (§ 10). Thus 




(116) 


Furthermore, at the top and bottom of the swing the reac- 
tion due to acceleration just balances the elastic restoring 
force of the shaft — that is 


— coV = kr 
g 


In fact, this relation affords a method of directly calculating 
the frequency. 

The equilibrium at this instant is just the same, how- 
ever, as that which would exist all the time if the weight 
W were to whirl in a circle of radius r at the same frequency. 
In the first case, the motion is given by y = r sin coL In 
the second case, the circular motion is given by the simul- 
taneous equations 


y = r sin co/ 
X = r cos (ct 


(117) 


When a rotor is revolving at such a speed, called a 
critical speed, a whirl is readily built up which may assume 
danprous proportions if the centrifugal force is aided by 
a slight unbalance. If it is speeded up above the critical, 

63 



64 


CRITICAL SPEEDS OF ROTORS 


a cliange of phase takes place which is discussed in § 47, so 
that the rotor becomes stable again. 

The frequency of whirl at which outward centrifugal 
force and restoring force e.\actly balance has been just 
shown to be identical to that for a linear vibration in a case 
like that in Fig. 41 for a perfectly balanced rotor. 

This same reasoning may be 
extended to the general case of 
a rotor of any mass distribu- 
Fig. 41 . -Shaft and FI>'^vheel. tion and Stiffness. If its vibra- 
tion frequencj’^ in a plane is 
found, its whirling speed or frequency is the same in all 
ordinary cases. It follows that the critical speeds of a 
uniform shaft can be found by the. same analysis as that 
used for flexural vibrations of bars in § 38. 

When once the transverse vibration frequencies of a 
rotor are found, however it be loaded, its critical speeds at 
which whirling may take place are approxi- 
mately known. 

46. Gyroscopic Action. — One effect, 
which enters for revolving rotors but is 
absent in a plane vibration, causes some 
difference between the linear and whirling 
frequencies in certain cases. This is a 
g 3 'roscopic effect due to the angular mo- 
mentum of the rotor. 

This effect does not come in near the 
center of span of the rotor, where the shaft 
is alwaj'S parallel to the bearing center line 
regardless of how much it deflects, as at 
A, Fig. 42. At B, things are different. 

The disk B has a so-called precessional 
upon it as the shaft whirls, which causes a restoring force 
additional to that due to the elastic stiffness, and thus 
produces a slightly higher frequency. Thus in a whirling 
at critical speed produced by unbalance the gyroscopic 
action tends to produce critical frequencies slightly higher 



Fig. 42. — Gyro- 
scopic Action of 
Flywliccl. 

motion forced 
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than frequencies due to linear vibrations. This effect is 
usually negligible and never larger, say, than 2 or 3 per 
cent. 

A simple experiment shows very well how the gyro- 
scopic force acts in this case. If a heavy top spinning at 
a high speed, at an angle of tilt such that it processes 
slowly owing to gravity, is grasped by the upper end of its 
axis and an attempt made to hasten its precession, it will 
show a strong tendency to assume an upright position. 
This is just what happens with the disk B. When a pre- 
cession is imposed upon it through the whirling motion of 
the shaft it tries to move inwards, thus producing a stiff- 
ening effect on the shaft. 

It is to be noted here that the opposite action takes 
place if the precession of a top is retarded or if it is made 
to precess in the opposite direction. It tends to fall over. 

On examining the direction of rotation of the top it will 
be found that forcing the precession in its direction of 
rotation causes the top to try to rise, and vice versa. 

In our problem the direction of whirl and the direction 
of rotation of the shaft are always the same, so the first 
effect is present. 

For a more complete mathematical treatment of this 
gyroscopic action see reference 8. 

47. Phenomena at Critical Speed in Detail. — From the 
foregoing reasoning the critical speeds of shafts may be 
found, but in order to build up a whirl a stimulus of some 
sort must be present, as in the case of the building up of a 
linear vibration. The usual stimulus responsible for 
building up a whirl is caused by unbalance, that is, the 
center of gravity of the rotor is not exactly on the axis of 
rotation. 

A detailed analysis will now be given of the simplest pos- 
sible case where a flywheel of mass m is carried at the 
middle of the span of a weightless shaft. In order to 
eliminate the consideration of gravity, the shaft will be 
taken in the vertical position. 
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Gravitation has no effect on critical speeds, however, 
and in the case of a horizontal rotor it cancels out in the 
analysis. 

Take a plane section through BC, Fig. 43, and take the 
XY plane in this plane. The origin 0 is taken on the 
center line betAveen the bearings to which S, the geometrical 
center of the flj^vheel, will return if the shaft has no 
elastic deflection. The elastic deflection is given by r, 
and a is the distance between the geometrical center of 
the fl>nvheel 5 and its center of gravity G. The distance 
a is actually very small, 2 or 3 mils, say, and equals zero 


B 



Fig. 43. — Unbalanced Flywheel on Revolving Shaft. 


for perfect balance. It is a constant for a given amount 
of unbalance as is assumed in this analysis. 

From Fig. 43 the acceleration of G is seen to be gov- 
erned by the following two independent equations, one 
for the X component of motion and one for the Y com- 
ponent. 

d~ r . , . . , dx 


m — [-v + a cos (cof -f a)l = — 6 -7 kx 

at- at 


VI 


df 


[y -f a sin (w/ -}- «)] = — — ky 

at 


(118) 


where is the friction or damping coefficient and k the 
elastic constant of the shaft. 
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These two equations are independent of each other, 
their solutions giving the X and the Y components of 
motion respectively. The first may be put in the form 

dH , ^ — inactP' cos (co/ + a) . (119) 

dfi dt 

Referring to Fig. 43 it is seen that the phase angle of r is 
taken as zero so that x has its maximum value when 
at = 0- 

The usual way of writing an equation of the type of 
(119) is to omit a, thus making the phase angle of the force 
maoP' zero instead of that of r. This causes r to lag behind 
mao? by the angle « as appears from the solution of (120) 
given by (123). There is a transient vibration term and 
a steady forced vibration term, but since the transient 
soon vanishes it is only the forced [vibration which is of 
interest. Writing the equation without a. it has the form 

m ^ ^ ^ ~ maoP' cos at . . (120) 


Solving in terms of mechanical impedance, we have 

ma<ir 

^ k — ina^ + jba 

where x and maoP are vectors. In like manner 


(121) 


y = 


jmaoP- 

k — ma^ + jba 


( 122 ) 


where the numerator contains a j factor because y must 
lead X by 90°, as can be seen from Fig. 43. Changing from 
the vector to the linear form of expression, as in § 28, 
these expressions become 

maa^ cos {at — a) 

V(^ — ma^Y "b 

maa^ sin (at ~ a) 

^ ~ V(k - ma^Y + 


. ( 123 ) 
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where 

bm 


These two expressions are 90° apart in both time phase 
and space phase, so that the motion is circular. 

As previously pointed out, the resonant frequency of 
the whirl is the same as if it were a linear vibration, with 
the stimulating force acting in one plane. 

48. Discussion and Interpretation of the Preceding 
Analysis. — Note that the preceding analysis assumes the 
frequency factor w to be a constant. The whirling is con- 
sidered as starting at a definite fixed to. This solution does 
not cover the cases for speed rising up through the crit- 
tical or falling off through it, solutions of which would be 
very desirable for a more complete study. The differential 
equations for these cases have not been completely solved 
(reference 9), whereas for « taken as a constant the solution 
is seen to be quite simple. 

The physical interpretation of the whirling motion of 
such a shaft at various frequencies is very interesting. 
Note that for a given value of w the angle a is fi.xed and 
OSG revolves about 0 as a unit. The term whirl as here 
used refers to a whirling motion of the shaft always coin- 
cident with the rotational speed co, as indeed it must be, 
because the vibration is forced by the unbalance and there- 
fore must be in period with it. There are other types of 
shaft whirling, where the whirling is always at resonant 
speed, regardless of the rotational speed of the shaft. 
These types of whirling are discussed in later sections 
(§§ 81 and S3). It is true, however, that whirling due to 
unbalance is severe only when the rotational frequency is 
in step with the resonant or natural shaft frequency, or 
ver}”^ nearly so. This is well shown in the analysis just 
given. At low speeds, a is small, r and a being almost in 
line and r being comparatively small. As resonance is 
approached, r rapidly increases, as does a, until at the crit- 
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ical speed is a maximum and a almost exactly 90°. As the 
speed rises still higher, r rapidly decreases and a. approaches 
180° so that G revolves inside instead of outside of S, and 
at very high speeds the system revolves almost exactly 
about G. At speeds high above the critical the tendency 
for the system to revolve about its center of gravity be- 
comes so powerful that the stiffness of the shaft is almost 
completely overcome, so that S is forced about G as a. 
center, just as a button spinning on a string will revolve 
about its center of gravity regardless of how the string 
passes through the button. When the shaft is running 
below its critical speed, centrifugal force which throws the 
center of mass outwards has the upper hand. When 
revolving above the critical speed the tendency for a 
revolving body to revolve about its center of gravity pre- 
vails. The critical speed is the transition point. 

Interesting experimental studies of the vibration phe- 
nomena at this transition point have been made which 
verify the analysis of the preceding section (see refer- 
ence 10). 

49. Dunkerley’s Formula for Critical Speed Frequency. 
— A very important early contribution to this subject of 
some forty years ago, which is a classic, was made by 
Dunkerley in England (reference 11). Experiments were 
made on model shafts carrying different numbers of fly- 
wheels and supported in different ways, which were checked 
against carefully worked-out analysis. The best-known 
result of this work is the following approximate formula 
for the critical speeds of shafts carrying a number of con- 
centrated masses or flywheels. 


2 W 2 


~ etc. . (124) 

Wr W32 


where co is the critical speed of the system ; w, that of the 
shaft by itself ; coi that of flywheel 1 alone on the shaft, the 
latter considered as weightless ; W 2 that of flywheel 2 alone, 
etc. Although not exact, this formula is very useful in cases 
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where there are two bearings as shown in Fig. 44. For 
three-bearing sets the use of this formula is to be avoided, 
as it may involve a large error. 

50. Energy Method — ^Rayleigh’s Principle. — The usual 
waj' of calculating the frequency of whirl of a rotor is to 
employ the principle of conser\'ation of energy. A plane 
vibration of some convenient amplitude is assumed, the 
maximum potential and kinetic energies estimated, and 
since these two must be equal to each other they are 
equated and the frequency constant w is solved for. 

To estimate these energies some form of deflection 
curve of the rotor must be assumed in which it vibrates, 
and the amount of the potential and kinetic energies 
will depend upon the deflection curv^e. It is at this 
point that an important general principle in vibration 

first proved by Lord 
Rayleigh will be em- 
ployed. Rayleigh’s 
principle may be stat- 

Fig. 44. — Shaft with Several Fi>-whccls. ed as follows: III any 

vibrating elastic system 
the displacement amplitudes are such that 'when the maxi- 
mum potential and kinetic energies are equated to each other 
the lowest or gravest possible frequency is obtained. If 
any other deflection curs^e than the correct one is taken, 
the frequency will come out a little too high. The system 
automatically deflects in such a way that its frequency is 
a minimum. 

It can be seen that this principle must be true, since a 
deviation from the actual deflection curv'^e of a vibrating 
system or shaft can be produced only by the presence of 
some additional elastic stiflfness to maintain the assumed 
deviation. Removal of material to produce the deviation 
is not permissible because this implies removal of part of 
the original system assumed. Any additional stiffness 
tlius produced, or restraint as it is called, increases the 
elastic potential energy and thus raises the vibration fre- 
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nuency. For a rigorous analytical proof of this principle 
see Rayleigh’s “Theory of Sound,” Vol. 1 (reference 12). 

Rayleigh also pointed out that if the static deflection 
curve of a loaded bar is used as an approximation instead 
of the actual vibration cur\'e, in calculating the energies, 
the error in frequency will be small. This is true because 
considerable change in load distribution is necessary to 
produce appreciable change in the deflection cur\’-e, so 
that the cur\'^e for vibration load distribution in a bar or 
shaft is nearly the same as the cur\^e for static load distribu- 
tion. It will be recalled from § 38 that the loading of a 
uniform bar in producing a given deflection cur%^e depends 
on the fourth derivative of that deflection curve, which 
explains the similarity of deflection curves for two quite 
different loadings of such a bar. 

51. Practical Use of the Energy Method.—In nearly 
all practical work the static deflection curve of a rotor is 
used in calculating its critical speed. For instance, for 
the shaft shown in Fig. 44 carrying the four flyAvheels, if 
the weight of the shaft Is neglected, the potential en&vgy 
of bending due to gravitation equals 

^W\yi -{- ^W 2 y 2 T etc (1-25) 

where yi, y 2 , etc., are the deflections due to gravity at 
Wi, W 2 , etc. The maximum kinetic energy of a vibration 
of this same amplitude equals 

co^yri-h^ — w^yo^ + etc. . . (126) 

2 g 2 g 

since 

. dy 

y = yi sm oil; -yr ~ wyi cos w/, etc. 
at 

Equating the potential to the kinetic energy and solving 
for / = w/27r 

f = / Wiyi -f K%2 + Wzyz W^y^ 


(127) 
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In the case of rotors with distributed loads the load 
deflection cup>-e is found by graphical statics, the shaft 
being arbitrarily divided into a number of sections for this 
purpose. The weights of these sections and their corre- 
sponding deflections are used as just shown, giving the 
general formula 


= /-IKv 


(12S) 


where the summation sign 2 is used for brevity. 

When required, a more accurate deflection curv^e can 
be found by using the static load-deflection curv'e to cal- 
culate a second one where the acceleration reaction of each 
section (Wui-y/g) is used in place of its weight. Anj' con- 
venient value of u is assumed, and the values of y are 
taken from the first curv^e. By Rayleigh’s principle the 
frequency thus found will be slightly lower than that which 
would be obtained by using the original static deflection 
cuiA'C to estimate the energ>'. 

In three bearing sets this process may fail in that the 
successive approximations give new frequencies which 
diverge from the true or gravest frequency rather than 
approach it (see § 55). 

52. Elasticity of Bearings and Critical Speeds. — Thus 
far, perfectly rigid bearing supports have been assumed 
in critical speed calculations. Actually, this ma)' be far 
from the truth, especially for large, heavy rotors. For 
every rotor there is a certain degree of elasticity in the 
bearings thcmseh'es, which lowers the critical speed below 
that for rigid bearings. This lowering of critical speeds is 
frequently as much as 25 per cent and may in extreme 
cases be over 50 per cent; tliat is, the critical speed may 
be reduced to less than one half that for perfectlj'^ rigid bear- 
ings. Furthermore, bearing pedestals will normally deflect 
sideways more easily than up and down, with the result 
that a rotor may have two distinct critical speeds, one of 
them a sidewaj's vibration and the other up and down, the 
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sideways vibration being the lower of the two, perhaps 
10 or 15 per cent lower than the other. The actual motion 
in such cases is not linear but elliptical in character, with 
the major axis several times the minor. ^ 

In critical speed calculations the effect of the yielding 
of the bearings may be taken account of by estimating 
their elastic deflections a and b as produced by the bearing 
reactions during vibration, and measuring the displace- 



Fig. 45. — Effect of Elasticity in Bearings. 


ments yi, ya, etc., from the line A'B' instead of AB as 
shown in Fig. 45. If the line of the shaft axis is e.'^tended 
from the bearings A and B, it will intersect A'B' extended 
at Ni and N 2 , which points may be considered the nodes 
of this type of vibration. The more the bearings yield, 
the farther apart are Ni and N 2 t and thus the lower the 
frequency. 


W, 


kj 




Q. 




Fig. 46. — Rotor with Elastic Base Supports. 


53. Effect of Elasticity and Mass in the Base. — If the 
rotor is mounted on a platform or structure which is itself 
not sufficiently rigid, critical speeds will enter which are 
characteristic of a compound vibrating system. This 
condition is shown in idealized form in Fig. 46. The 
weights and elasticities of the system are marked to corre- 
spond with that analyzed in § 33, the stimulus, in this 
case, coming from unbalance. As in the former instance. 
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the system has two natural frequencies: the lower where 
Wi and TT'o vibrate in phase; and the higher where they 
vibrate against each other, 180° out of phase. The latter 
is likely to produce the sharpest shaft deflection, although 
the possibility of its presence is often overlooked in design 
as is also the compound character of the vibrating system. 
In studying the effect of mass of and elasticity in the 
base on these two critical speeds it is instructive to plot 
cuiA-es showing the variation of the critical speeds with 
variation of the ratio W 2 /W 1 and ^ 2 /^ 1 . 

Thus far, only the up-and-down criticals have been 
considered. If the base is practically rigid sideways, only 
the regular sideways shaft critical speed is present in addi- 
tion, making three critical speeds in all. 

54, Upper Critical Speeds. — Every shaft has a series 
of critical speeds, of which the first or fundamental 
is much the most important. A continuous elastic rotor 
like a shaft has an indefinitely large number of higher 
critical speeds according to ideal theory, but in practice 
only the first two or three are of importance. Those 
higher than this are relatively difficult to excite, and are 
almost always well above operating speed. Upper critical 
speeds have not been mentioned thus far except for the 
compound system of the previous section. 

The characteristic deflections and frequencies for the 
more important critical speeds of a uniform shaft sup- 
ported in two and in three equally spaced bearings will 
now be compared. 

Figure 47 shows the deflection curves of a uniform 
shaft supported between two bearings for the first and 
second critical speeds. These deflection curves and their 
frequencies are found by the same analysis as that given 
for flexural vibrations of bars presented in § 39. The 
cur\'es in this case are of sinusoidal form, the first critical 
having one loop with nodes at the bearings, the second two 
loops with a node half way between the bearing nodes, the 
third with two nodes equally spaced between the bearings. 
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etc. The frequency is in each case proportional to the 
square of the number of loops, that is, the second critical 
is 4 times the first, the third 9 times the first, etc. 

For a uniform shaft supported in three equally spaced 
bearings the first three critical speed deflection cur\^es are 
as shown in Fig. 48. 

In this case, the first critical speed is the same as for 



Fig. 47. — Critical Speed Vibrations of Shaft. 


a two-bearing set of half the length of the three-bearing set, 
as though the shaft were cut at the middle bearing. This 
must be true because the middle bearing is at a node where 
both the flexure and the reaction due to vibration are zero. 

The second critical is a direct result of the presence of 
the middle bearing which supports a vibration reaction 
and at which there is a point of max:imum flexure. 


First critical 




Seconal critical 
16 times first 



Third critical 
4 times first 


Fig. 48. — Critical Speed Vibrations with Tliree Bearings. 


^ The third critical corresponds exactly to the second 
critical for a single span just as the first corresponded to 
the first for a single span, as previously noted. Like the 
first critical speed, the frequency of the third would not be 
altered if the shaft were cut at the middle bearing. Thus 
the frequency of the third critical is 4 times that of the 
first, but that of the second about 1.6 times the first. 
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The frequency of the second critical for such a shaft in 
three bearings can also be calculated from the anal 3 ?sis of 
§ 39, using the conditions at the middle bearing in addition 
to the ends in eliminating the constants in equation (91). 
For instance, referring to Fig. 48, at the middle bear- 
ing u = 0 and the slope dyjdx — 0. The slope, however, 
would not be zero at the middle bearing if the spans were 
not of equal length, so the solution then would not be 
so simple. 

55. Practical Application. — The ideal cases presented 
in the previous section are useful for illustration but are 
seldom realized in practice for several reasons. The mass 
distribution of a rotor is not uniform, but is concentrated 
along the middle part of the span. For two-bearing sets 
and also for three-bearing sets this has the effect of bring- 
ing the criticals previously discussed closer together. For 
relatively heavy flywheels and light shafts they may be 
quite close together. -Three-bearing sets have the com- 
plication that unbalance causing vibration may be present 
in both spans and not in the same plane. For this reason 
three-bearing sets are often troublesome to balance. 
With four bearings it is sometimes possible to shift the 
rotor of one span with respect to the one in the other by 
uncoupling and coupling again at a diflFerent angle so as 
to bring the unbalance into the same plane, and thereby 
facilitate balance. The two middle bearings of a four- 
bearing set are usually so close together that such a set 
acts like a three-bearing set as regards critical speeds. 

In such cases, the energy method is used to find the 
critical speeds as presented in §51, but for three-bearing 
sets the true deflection cur\'^e is often difficult to establish 
even approximately, and errors must be looked out for. 
Nearly all the vibration amplitude may be present in 
one span if the span is long or carries a particularly heaA^ 
rotor. Figure 49 shows a typical case. 

56. Reciprocal Theorem . — h very interesting general 
theorem in elasticity first published in simple form by 
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Maxwell in the Philosophical Magazine about seventy 
years ago (reference 13) will now be introduced because, 
in addition to its general interest, it is sometimes used in cal- 
culating the critical speeds of three-bearing sets to find 
the reaction on the middle bearing. 

If two separate points on any elastic body be selected 
and a given weight placed at the first point, the deflection 
of the second point thus produced is the same as the deflec- 



Fig. 49. — Unsymmetrical Vibration at Critical Speed. 


tion produced at the first point on placing the weight on 
the second point. For instance, a force F acting at P\ 
on the cantilever shown in Fig. 50 produces the same down- 
ward deflection at P 2 as that produced at Pi when the 
force is removed from Pi and acts at Pg. 

The proof of the reciprocal theorem for the particular 
case of this cantilever of Fig. 50 will now be given. 

Let yi and yo be the deflections produced at Pi and Po 
when P is applied alone. 



Fig. 50. — Reciprocal Theorem. 


^ Let y-i! and y^ be the deflections produced at these 
points when F' is acting alone. 

The potential energy of bending produced by first 
applying F by itself and then F' in addition evidently 
equals that produced by applying the weights in the 
reverse order, but at the same points as before, because 
the final deflection in each case is produced by the same 
two forces acting in the same positions. 
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Expressing the energy' mathematically 

F alone F' In addition F* alone F In addition 

at Pi at Pi at P* at /*i 

S-I + hF'y'. + Fy'i = hfy2 + hFyi + Fyo (129) 

Note that in applying the second force, in each case 
the first force already acting in full does additional work. 
Cancelling off identical terms from each side of the equa- 
tion gives 

Fy\ = F'yo (130) 

This is a mathematical statement of the reciprocal theorem 
for the case of the cantilever of Fig. 50. 

If F and F' are equal, as is the case when the same 
weight is used in each position as assumed in the first 
illustration, 

y’l = y2 

which shows that if the same weight is used in each posi- 
tion the two alternate deflections are exactly equal. 

The proof given is seen to assume that the deflections 
produced on adding a force are the same, whether or not 
initial deflections, due to some other force, already exist. 
This illustrates another fundamental law of elasticity 
which follows from Hooke’s Law. It is called the principle 
of superposition of strains and holds to a close approxima- 
tion provided the strains are not so great as to materially 
change the shape of the body, such as an excessive down- 
ward bending of the cantilever or the large deformations 
of a rubber body. The whole mathematical theory of 
elasticity assumes the truth of the principle of superposi- . 
tion of elastic strains, whereby a given loading produces 
its own particular strains on top of strains already present 
from previous loadings as though these initial strains did 
not exist. 

57. Application to Three-bearing Sets. — ^As an exam- 
ple, assume a three-bearing set loaded as shown in Fig. 51. 
The problem is to find the middle bearing reaction Rb. 
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First, calculate the deflection j'b for a load LatB with 
the bearing removed. Then, from Hooke’s Law that strain 
is proportional to load, if Rb produces the deflection yB 
(with the bearing B removed), it follows that 


Rb _ Jn 

T~Yb 


(131) 


In this relation, if yB is known, Rb can be calculated, 
since L and y ' b are known. Now y'B can be found from 
Maxwell’s theorem as follows. Again, with bearing B 
removed, if L produces deflection yi at 1, 3*2 at 2, etc., 
unit load at B produces deflections yi/L, Ji/B, etc., respec- 
tively. 



Fig. 51. — Forces on Shaft of Tliree-bearing Rotor. 


By Maxwell’s theorem, unit load at 1 produces deflec- 
tion ji/L at B; at 2 it produces y^jL at B, etc. It fol- 
lows that load Fx at 1 produces FiyxIL at B, and load Fz 
at 2 produces deflection FzyzjL at B, etc. Thus the total 
deflection at B = {'SiFy)fL. But this equals yB, since 
the bearing load just compensates this deflection, that is 

y^= (132) 

Substituting (132) in (131), 

Bb = ^ (133) 

y B 

where F has the values Fi, Fz, etc., and y the values 
y 2 , etc., produced by L at these points, and y'n is the 
deflection L produces at B. The problem is thus solved 
by calculating the deflection curve between A and B for 
a single arbitrary load at B and by appljdng Maxwell’s 
theorem. 
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TURBINE WHEEL VIBRATIONS 

58. Flexural Waves in Turbine Wheels. — During the 
few years before 1920, enormous development took place 
in steam turbines, which has indeed continued since then, 
but about that time a number of disk wheel failures took 
place Avhich were not understood. An investigation co^"- 
ering a period of two years proved beyond doubt that these 
failures Avere due to fatigue cracks produced by a particular 
type of axial resonant vibration. When the presence of 
this resonance Avas once demonstrated by tests, its cause 
Avas easy to picture, as is often true in such cases. 

It is AA'ell knoAA’-n that sound Avaves traA'^el AA'ith a definite 
A'^elocitj' characteristic of the density and elasticity of the 
medium Avhich transmits them. In like manner, the edge 
or periphery' of an elastic disk Avill transmit a train of 
flexural AvaA'es around it Avith a characteristic speed deter- 
mined b)’’ the mass per unit area, and the stiffness of the 
disk. They differ from sound AvaA^es in that they invoh^e 
a bending or flexure of the elastic material instead of dirept 
compression and rarefaction. Since the boundary of the 
disk has a finite circumferential length the AA'aA'^e lengths 
must haA^e certain fixed A'alues. For instance, the disk 
edge may carry a continuous train of 2, 4, 6, or S, etc., 
AA-aA-es Avhose AvaA^e length is v, j, J, I, etc., of the circum- 
ferential length of the disk. Furthermore, instead of all 
AA'aA’e trains moAung Avith the same characteristic velocity, 
as in sound, each of these AA'aA'-e trains has a different 
A’elocity, the longer AvaA^e lengths being transmitted AA'ith 
considerabl}'- higher speeds than the shorter. 

59. Critical Speeds of Turbine Wheels. — It Avas found 
that Avhen the forw'ard rotational speed of a turbine Avheel 
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was exactly the same as the characteristic speed of one of 
the wave trains a dangerous condition existed because the 
wave train was easily built up by a comparatively small 
excess of axial pressure on some part of the wheel. The 
wave train always moved along the disk edge against the 
direction of rotation of the wheel in such a way that it 
remained stationary in space owing to the fact that its 
backward speed was the same as the forward rotational 
speed of the Avheel. A little consideration will show that 
a single pressure spot on the Avheel will build up such a 
wave train if it acts in the right phase relation with the 
wave. If in Fig. 52 the bottom arrow represents the wheel 
motion and the top the wave motion of equal velocity, the 
waA^e shoAATi aauII remain stationary on the diagram as the 
Avheel edge moA'^es to the right. A fixed force at P acting 



Wheel velocity 

Fig. 52. — Wave in Turbine Wheel. 


betAveen the crest and trough of the AAmA'^e AAdll continually 
supply energy to the AAmA^e train. The rate of energy 
supply equals the Avork per second AA^hich equals F times the 
transverse A^elocity component in the direction of F. Let 
the Avave equation be 

y = 3’o sin cot = yo sin l-n-ft 

Avhere/ = frequency. 

If X = the Avave length, / = F/X, so that 

. 2TrF/ 

y = yo sm -y- ( 134 ) 

This equation gives the transverse vibratory motion 
of a particle of the AA^heel as the particle moves to the right 
against the standing Avave. The transverse velocity of the 
particle = dy/dt = 27ryoF/X cos 27rF^/X. As the particle 
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passes through the fixed point P it has its maximum trans- 
verse velocity = 27ryoTV^- 

The force P at P acts continually on a succession of 
such particles, all of which have the -transverse velocity 
liryoV/'K. Thus 

Work per second = ° . . . (135) 

X 

This rate of energy supply is easily found for a given 
velocity of motion V, wave length X, and force F, for any 
amplitude yo. 



0.10 0.20 0.30 0.40 

Total amplitude, inches 

Fig. 53. — Energy Supply and Dissipation in Vibrating Turbine Wheel. 


Calculation shows that, if F is even a few pounds in 
magnitude, long waves on a large, high-speed wheel are 
easily built up. 

The curves of Fig. 53 show the energy supply curve 
and the energy dissipation curve as dependent on ampli- 
tude for a typical case, in watts. The wave amplitude 
builds up until the vibration dissipating power of the disk 
equals the energy input, where the curves intersect. 

Note that the ability of such a large steel wheel to dis- 
sipate energy'- is comparatively small, and that only a few 
watts of energy' input may produce a dangerous vibration. 
The dissipation curve shown was obtained from a direct 
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measurement. The input curve is based upon formula 
(135), assuming reasonable values of F and V. Note also 
that the input curve is linear, but the dissipation curve 
shows an increasing rate of dissipation with amplitude. 
Thus at very small amplitudes the probability of a wave 
building up is large, since the input curve is well above the 
dissipation curve in this region. 

60. Prevention of Turbine Wheel Vibration. — Such a 
critical speed vibration of a turbine wheel is a serious danger 
even if the wave amplitude be comparatively small, because 
of the possibility of a fatigue crack developing in the wheel 
from the periodic stresses present, often resulting in failure 
if not discovered in time. It is now known that many 
steels are much more susceptible to fatigue failures in the 
presence of steam at high temperatures, which must be 
taken account of in turbine wheel design. The allowable 
stress under these conditions is many times less than that 
permissible in a structure under static load. 

Just as in shaft and rotor design, the means employed 
to prevent turbine wheel vibration is to make sure that the 
critical speeds are different from the running speed. The 
problem is much more difficult in the case of turbine wheels, 
however, because the dangerous critical speeds are numer- 
ous and hard to avoid. In very large turbine wheels it is 
impracticable to design them with all critical speeds above 
the running speed because of their large size and mass, so 
the critical speeds are sometimes spaced on both sides of 
the running speed. 

The flexural wave trains responsible for the critical 
speeds of course traverse the wheel edge in its completely 
bucketed condition, the wheel with its buckets vibrating 
as a unit. The most dangerous critical speeds are those 
involving two and three waves or the so-called four- and 
six-node criticals, because the waves are long and extend 
radially inwards below the bucket zone, so that the maxi- 
mum flexures are in the wheel itself, involving the danger 
of wheel failures. Where the number of waves in the train 
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is large, only the buckets are apt to be involved, the great- 
est flexures occurring in the bucket zone, with danger of 
bucket failures, hlany obscure cases of bucket failures 
have been traced to a train of many short wa^^es carried 
along in the bucket zone well outside of the wheel proper. 
In such cases there may be as manj'' as eight or ten waves 
in the entire circumferential train. 

Although the wave velocities are susceptible to accurate 
calculation from theoiy', it has been found that, for safety, 
tests and frequency' measurements must be made on 
revolving wheels on a test shaft in a steam box. Theory 
does not take account of conditions of initial stresses, 
which are present to a greater or less extent in the best 
turbine wheels. It is easily demonstrated bj'’ test that a 
thin disk with initial hoop tension or compression along its 
edge will vibrate at radically different frequencies from 
those of an ideal stress-free disk. 

The effect of centrifugal force is to raise the flexural 
stiffness of the disk and increase the wave speed, but this 
effect is appreciable only for the longer waves, and further- 
more it is easily calculated. See reference 14. 
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61. Unbalance in a Rigid Rotor. — If the axis of rota- 
tion of a rigid rotor passes exactly through a principal axis 
of inertia (§ 42), its rotation will be perfectly smooth; but 
if the axis is slightly displaced from this ideal line, periodic 
bearing reactions in step with its rotational frequency will 
result. Figure 54 shows a cylindrical rotor whose prin- 
cipal axis about which rotation should take place for perfect 
balance, as it is called, is supposed to exactly coincide with 
the axis of the cylinder. Unbalance is assumed, however, in 
that the axis of rotation is displaced from this ideal line 



Fig. 54. — Example of Unbalanced Rotor (Exaggerated). 


as shown exaggerated in the figure. This illustrates the 
most general case of unbalance, where the rotational axis 
in its failure to coincide with the principal axis of the 
rotating mass does not even pass through its center of 
gravity. Practically all manufactured rotors, large and 
small, show more or less of this general type of unbalance, 
a very small amount of which causes vibration, , especially 
in high-speed rotors. Regardless of precision in manu- 
facture, slight irregularities in mass distribution are found 
to be present, sufficient to cause perceptible vibration, even 
though the shaft passes exactly through the geometrical 
axis of the rotor. 
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That part of the unbalance due to the axis of rotation 
not passing through the center of gravity of the rotor 
is termed sialic unbalance, because it can be detected by 
a delicate static test. If the rotor is placed with its bear- 
ings on smooth, level ways it will roll under the influence 
of gravity until its center of gravity is in the downward 
position. That part of the unbalance which remains due 
to a slight lack of alignment of the principal axis and tlie 
rotational axis is termed dy7iamic unbalance, because it 
can be detected . only by a rotational test bringing into 
play the dynamic forces of rotation, or the centrifugal 
forces. 

62. Analysis of Unbalance of a Rigid Rotor. — It might 
be gathered from the preceding section that two steps are 



Fig. 55. — Complete Balance Established by Two Weights. 


necessary in the balance of a rotor, one for the static and 
one for the dynamic unbalance. This is not so, however, 
as all unbalance present must show up in a rotational test. 
It will now be shown that only two balance weights are 
necessary to correct the unbalance of a rotor in the most 
general case. 

The rigid rotor in Fig. 55 illustrates such a general case 
of unbalance. The center of gravity is displaced from the 
axis of rotation by the central weight Wo opposite the 
center of gravity while at the same time a couple is pro- 
duced by the pair of weights Wi and W 2 in a different lon- 
gitudinal plane from Wo- It will be recalled that a couple 
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tends to make a free body turn about its center of gra^^ty 
regardless of where the couple is applied (§ 40). Thus an 
unbalanced couple in the axial plane of a rotating body, 
regardless of where the couple acts, may be represented by 
TTi and W 2 as showm in the sketch. The force vectors 
represent the reactions produced by the weights during 
rotation. The force vector Fq is replaced by its equivalent 
of two component vectors F' and F” in the transverse 
planes of Fi and Fo respectively. The vectors F' and F^' 
are the equivalent of Fo if they have equal and opposite 
moments about the center of gravity of the rotor and at 
the same time the sum of their magnitudes equals that of 
Fo. If the resultant vectors Ri and Ro are each exactly 
counteracted by the balance weights Bi and B 2 , complete 
balance will be established. Any case of an unbalanced 
rigid rotor may be represented in this way. 

Note that TFo produces the static unbalance, and TT'^i 
and W 2 produce the centrifugal couple which is the cause 
of dynamic unbalance and can be detected only by a rota- 
tional test. Note also that a rotational test is sufficient 
to detect all unbalance in the rotor. 

63. Methods of Balance. — Modem types of balancing 
machines employ what may be called the double fulcrum 
method. The rotor is -elastically supported in a light, 
stiff frame which is first pivoted in the plane of Bi while 
weight B 2 is applied, and then pivoted in any other plane 
while weight Bi is applied (§ 62). By this method the 
addition of Bi has no influence on the compensating action 
of Bo because when the frame is pivoted in the plane of 
Bi the moment of that weight about the pivot is zero. 
\^ffien the moments about two separate planes are at the 
same time zero, no moments tending to vibrate the rotor 
can be present, and complete balance is thus established. 
Figure 56 shows the plan of a simple outfit for balancing 
small rotors, using the traverse method. The rotor is 
placed in the frame F, which for small rotors may be of 
wood, as such a frame can be made light and stiff. The 
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rotor R rests in half bearings in the frame and is driven by 
the motor M through the belt L whose tension is adjusted 
by the screw 5 by which the motor on its platform P can be 
slid to and from the rotor. The rotor is first pi^•oted in the 
plane C, say, by the adjustable pivots, those in the plane D 
being withdrawn. A spring, not shown, is placed under the 
frame at the bearing B, of sufficient stiffness to cause reso- 
nance at about 700 or 800 r.p.m. The rotor is driven up 
through this resonant frequency, and the deflection indica- 
tor E is adjusted with its button resting on the top of the 
frame so that it registers the vibration present. The arc 
through which the needle oscillates is an appro.vimatc 



measure of the unbalanced moment about the fulcrum in 
the plane C. The indicator should read down to mils 
(thousandths of an inch) at least. A standard Starrett 
indicator sers'^es well. 

To find the amount and position of the weight the 
\nbration amplitude is first read and then a trial weight 
of small size is successively placed in four positions 90° 
apart in the plane D, the amplitude of vibration shown by 
the indicator being read in each case. A cun,'e of rise and 
fall of amplitude as dependent on the angular position of 
the trial weight is then plotted, from which the amount 
and position of the weight can be estimated. An inspec- 
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tion of the curve shows in what position the weight sub- 
tracts most from the unbalance and Avhere it adds most 
to it, the former position being where the balance weight 
belongs. Its amount can be estimated by noting the 
change of vibration amplitude produced by the weight in 
different positions and comparing with the original vibra- 
tion amplitude present, assuming that the vibration ampli- 
tude is a direct measure of the weight required. 

A piece of wax may be used as the trial weight for 
small rotors. Besides having a light frame to insure large 
amplitudes, the outfit must not vibrate as a whole because 
of insecure foundation, nor must it be placed where floor 
vibrations are transmitted to it, otherwise inconsistent 
readings will be obtained on the indicators. A good way 
to isolate the appara- 
tus from such vibra- 
tions and prevent 
detrimental shaking 
is to place the outfit 
on a massive iron 
plate which rests on 
toy rubber balls or Pig. 57, — Rotor in Balancing Frame, 
on an inflated inner 

tube of an automobile tire. Also, the driving motor must 
be well balanced, and a light belt used to give a constant 
tension, which must be sideways, that is, its pull should be 
perpendicular to the oscillatory motion of the frame, nor 
should the pull be too great. 

To balance the other end of the rotor, the process is 
repeated, using the other fulcrum. The weights thus deter- 
mined should produce complete balance so that no oscil- 
latory motion is shown on the indicators. When the 
process is thoroughly understood a little practice Avill 
insure great proficiency on the part of the operator, both 
as to speed and as to accuracy. 

^ 64. Analysis of Frame Vibration. — In Fig. 57 the frame 
F is shown pivoted at P and carrying a rotor. Let us find 
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the vibration produced by an unbalance represented by 
the weight Ws- Let the moment of inertia of the frame 
and rotor about P = I. Torque amplitude about P due 
to Ws is equal to 


L = 


WBU-rh 

g 


(136) 


U k — elastic constant of the spring S, the torque per 
radian of arc about P — L\ = kl<t, since h is the lever 
arm, and for a radian of motion the force on the spring 
would be kh', and Li = force times lever arm. Thus 


U = hh Xh = kh- .... (137) 


The amplitude of the angular oscillations about P is 
found by the mechanical impedance method in a way per- 
fectly analogous to that used in the case of a linear vibra- 
tion where F is replaced by L, in by I, k by Li, and the 
displacement y by the angle 6. Thus 


6 X impedance = B 




6 = 




g 



(138) 


If w is large so that is relatively a small term, it 
may be dropped so that 


6 = 


- Wsrh 
I 


(139) 


where the minus sign shows B and the exciting forces are 
in phase opposition, and may be dropped in considering 
magnitudes only. 

It is seen from this e.xpression that the angle B through 
which the frame oscillates is proportional to the amount 
of unbalance present, as assumed in § 63. Furthermore, 
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the amplitude angle is inversely proportional to the mo- 
ment of inertia of the system about P so that to make the 
indicator readings as large as possible in order to detect 
small amounts of unbalance the frame must have the 
smallest possible inertia consistent with stiffness. Rea- 
sonable stiffness of frame is essential, since deflections in 
the frame disturb the indicator readings. 

Many of the older commercial balancing machines have 



Fig. 58. — Thearle Balancing Machine. 


frames far too heavy for accurate work, because of tlieir 
large inertia compared with that of the rotor. 

65. Thearle’s Automatic Balancing Machine . — h very 
interesting and efficient type of balancing machine has been 
worked out by E. L. Thearle, of the General Electric Com- 
pany, the theory and operation of which are described in a 
paper presented before the Applied Mechanics Division 
of the American Society of Mechanical Engineers in June, 
1931 (reference 15). A brief description of it will be pre- 
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sented in this section. The double fulcrum method is 
used as described in § 63. Figure 58 shows a photograph 
of the machine. Balance is accomplished by two weights 
just as in the apparatus shown in Fig. 56, but the amount 
and position of each balance weight is found from a single 
observation of a special automatic balancing head {d, Fig. 
58, and Fig. 59). 

Referring to Fig. 58, the rotor is seen to be mounted 
in a frame which is made up of a pair of light steel tubes 
with adjustable aluminum cross- members in which half 
bearings are mounted. The whole frame is carried on a 
four-point spring suspension supported by the curved 
uprights shown. A pair of adjustable fulcrums are pro- 
vided, consisting each of a pair of phosphor bronze springs 
held at their lower ends in a pair of vises, one pair for each 
supporting plane, which can be clamped or released at will. 
Three of these vises, and also their clamping rods with 
handles at the ends, appear in the figure. Either pair of 
fulcrums is thus readily put into action or released by a 
twist of the proper pair of handles. The position of these 
fulcrums on the frame is also adjustable, depending on the 
size of the rotor and the required balance weight planes. 

The automatic balancing head is shown in Fig. 59. 
The collar A is clamped over the end of the motor shaft by 
giving it a half turn, this fixing the head in position. The 
rotor is then driven up through the critical speed of the 
balancing frame so that it seeks its own axis of revolution 
(§43). 

The balancing head is then whirled bodily by the shaft 
in such a way that, when the steel balls C are set free in 
their race by a pressure on G, releasing the clamp D, they 
immediately roll to a position in the race opposite to the 
hea^Y spot on the rotor, thus cancelling the effect of unbal- 
ance, and reducing the vibration to zero. If they roll too 
far the rotor becomes unbalanced the other way, inducing, 
in turn, a whirling vibration of the head such that the steel 
balls wll tend to move back part way, etc., until they 
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reach a position of stable equilibrium when perfect balance 
is attained such that the head spins exactly about its axis. 

Note that operation above the critical speed of the 
frame is necessar^q otherwise the rotor vibrates with its 
hea\'y side out, so that when the balls are released they 





(c) 

Fig. 59. — Automatic Balancing Head. 


assume a position which makes the unbalance worse instead 
of improving it. 

In counteracting the unbalance, the steel balls assume 
positions such that the vector sum of the centrifugal forces 
produced by them just counteracts the unbalance, being 
nearer together or farther apart according as the unbalance 
is large or small. With no unbalance they take positions 
opposite each other. If the unbalance is excessive they 
fail to eliminate unbalance even if they come together. 
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Nevertheless, the plane of unbalance is indicated, so that 
a trial balance weight can be attached which so reduces 
the unbalance that a second trial will give complete bal- 
ance. From the positions of the steel balls the amount of 
unbalance is determined in both phase and magnitude, and 
a correction weight can be easily calculated and placed on 
the rotor. 

Balance accomplished by machines of this type is rapid 
and accurate, as has been proved by their use in production. 

For very small or very large rotors, other methods of 
balance are desirable. For small rotors the manufacture 
of the head becomes a watchmaker’s job; furthermore, its 
mass becomes so great compared with that of the rotor 
that sensitivity is reduced. 

For large rotors of over a ton in weight other methods 
such as that of the Akimoff type of machine require a less 
expensive outfit and one which is sufficiently rapid for 
massive rotors which are not balanced in quantity. 



CHAPTER XI 


PREVENTION OF NOISE AND VIBRATION THROUGH 
ELASTIC SUSPENSION; THEORY OF VEHICLE 
SUSPENSION 

66. Noise Elimination Classification. — An important 
field under the subject of vibration prevention is the elimi- 
nation of objectionable noises and hums which have their 
origin in vibrating bodies. Sound itself is a type of vibra- 
tion, and we might go far into the field of acoustics in noise 
elimination work. In this discussion it will be touched 
upon only as far as it is necessary' to understand the prin- 
ciples at the basis of noise and vibration prevention through 
elastic suspension. 

Noises from machinery', from the standpoint of elimina- 
tion, fall into t^vo classes: 

(1) Those which are radiated directly into the air 

from parts of the apparatus itself. 

(2) Those which do not radiate in this way, but 

are transmitted through the foundations to 
the building structure and appear with objec- 
tionable intensity often at remote parts of the 
building such as walls and panels which may 
resonate to the exciting frequency. 

Preventive methods may also be classified as follows: 

(1) Prevention at the source by study of the design. 

(2) Elimination through screening. 

(3) Elimination through elastic suspension. 

The first of these measures is always aimed at but, 
even if at all possible, may require expensive investiga- 
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tions and more knowledge of design than is at present 
possessed. 

Of the two remaining preventive measures it is evident 
that air-bome noises are not to be eliminated by elastic 
suspension, but that sound-absorption screens and filters 
must be resorted to. 

Elastic suspension is a specific for noises transmitted 
through foundations, and this is the phase of noise elimi- 
nation which will be treated in the few sections which follow. 

In general, air-borne noises are of comparatively high 
frequency, say 500 cycles per second and above, whereas 
foundation hums are seldom of frequency above 500 
cycles, and are usually considerably below this figure, and 
may consist of very low frequencies indeed such as those 
arising from unbalance of rotating bodies. 

67. Theory of Elastic Suspension. — The purpose of 
elastic suspension is to prevent the transmission of vibra- 
tions to the floor or foundation on which the vibrating body 
is placed. In order to specify the effectiveness of an 
elastic suspension its transmissibiliiy must be known. 
This is best defined as the ratio of the disturbance pro- 
duced in the foundation when the spring suspension is 
used to that when the vibrating body is rigidly attached, 
usually specified by the Greek letter e (epsilon). C. R. 
Soderberg, of the Westinghouse Electric and Manufactur- 
ing Company, was the first to correctly formulate the 
fundamentals of this subject (reference 16). For pur- 
poses of analysis, e is defined as the ratio of the vibration 
amplitude produced in the foundation when using the 
elastic suspension to that produced with a rigid connection. 
In many simple cases the foundation is relatively rigid 
and therefore its motion is comparatively small, so that 
the amplitude produced is proportional to the force of 
reaction. The value of e can then be found easily as the 
ratio of the force transmitted to the foundation through 
the elastic suspension to that transmitted through a rigid 
suspension. 
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Suppose that the mass m supported by the spring k is 
acted upon by a periodic force of amplitude F as shown in 
Fig. 60. Then using the method of mechanical impedance, 


yi{k — moP) — F 
F 

= T ^ 

k — mw~ 


and the force transmitted to the founda- 
tion through the spring is 


kyi = 


kF 

k — mur 


y, 


m 


77//P7///F7/F/7> 


Fig. 60. — Elastic 
Cushion for Vibrat- 
ing Mass. 


That transmitted with k replaced by a rigid support is 
simply F. Therefore 


kF 


e = 


k — mor 
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F = 


k — Vtur 
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/c“ - P 


(140) 



where fc is the free vibration 
frequency of the mass m on 
the spring k. 

If ///c , 


1 . 


€ = 


1 - r~ 


(141) 


Figure 61 shows how e varies 
with / as given by equation (141) 

Fig. 6i.-Vibration Transmissi- for a given resonant frequency 

bility Curve. € being always counted positive, 

since the negative sign which 
appears tvhen ///^ is greater than unity merely signifies a 
phase opposition between the two forces, which does not 
concern us. It is evident that, for any reduction in trans- 
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missibility at all over that produced by rigid suspension, € 
must be less than unity. Experience shows that, in ordinary' 
cases of sound insulation, for an elastic suspension to be really 
effective, e must be at least as small as ^ to From 
Fig. 61 it is seen that if the spring suspension is so stiff 
that it resonates at a frequency such that fc is higher than 
the exciting frequency so that r is less than unity, e is 
always greater than unity and the suspension is worse than 
useless. It becomes effective only when fc is to y of /, 
that is, when the spring suspension is comparatively limber. 
For example, if f/fc = 7, 

e = ^ ^ (dropping the negative sign) 

which gives good cushioning. 

For the case of a torsional vibration the same analysis 
applies, giving the same formula for transmissibility, the 
elastic support being of course designed to cushion a tor- 
sional vibration in such a case, the symbols of equation 
(140) referring to torsional frequencies. 

68. Special Case. — ^An elastic support may have more 
than one degree of freedom, each having its own natural 
frequency, fc. Each of these values of fc must be suffi- 
ciently below the exciting force / to insure sufficient cush- 
ioning. If one of these frequencies happens to be near / 
that particular mode of vibration is very likely to be picked 
up to some e.xtent even though the exciting force / acts in 
a different direction from the line of displacement of that 
vibration. For instance, the present General Electric 
refrigerator pump mechanism is mounted on a three-spring 
suspension as indicated in Fig. 3. This has four definite 
independent frequencies to which it readily responds: 
(1) torsional, (2) sideways tilt, (3) up-and-down bounce, 
(4) sideways shearing oscillation. The last is highest of all 
and nearest the pump frequency /. All values of fc, how- 
ever, are far below the 120 motor hum, which is thus 
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well isolated. The pump frequency comes through some 
as a slight quiver which is not objectionable because it is 

noiseless. ^ -j 

69. Effect of Damping in Spring Supports. — Consider 

a case like that of § 67 except that the spring support con- 
tains a small amount of damping.* Figure 62 illustrates 
this case where the damping factor h acts in parallel with 
the elastic constant as though it were 
contained in the spring material. From \f 

the rules of § 31, 


771 
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yi{k — -\-jbo:) = F 
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k — nioi- -|- jb 
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Fig. 62. — Elastic Sus- 
pension with Damping. 
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Thus, expressing e as a magnitude and dropping the phase 
angle. 
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since w? — 'k/vi. 
namely. 


6 


This can be put in terms of frequencies, 


// + 


EIL 

47r^7;7." 


(p-pr~ + b^ 


p 


4Trhn- 


(143) 


This expression reduces to (140) when i = 0 as of course 
it should. 


* See reference 4. 
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Note that damping slightly increases transmissibility. 
The increase is small, even when rubber is used as an 
elastic cushion, although rubber contains considerable 
damping. Thus, contrary to the belief of many engineers, 
damping in an elastic suspension may increase vibration 
transmission. Furthermore, springs do not damp a vibra- 
tion, but prevent vibration transmission. Damping im- 
plies energy dissipation in the form of heat, and springs 
do not dissipate energy. Dissipation is effective, however, 
in reducing the amplitude of a resonant vibration in which 
the amplitude builds up to a peak value. Vibration ordi- 
narily is transmitted to the base of a machine as a forced 
vibration. If the foundation were truly resonant, suc- 
cessful operation would be almost impossible. 

70. Effect of Foundation Char- 


:;F acteristics on Transmissibility. — Im- 

m; I y, portant cases arise where the mass 

and elasticity of the foundation itself 
- 4 must be considered. The mechanics 

I of this case is represented in Fig. 63. 

k. The transmissibility is defined as be- 

fore as the ratio of the amplitude of 
Fig. 63.— ElasticSuspen- vibration transmitted to the founda- 
sion on Elastic Floor, tion through the spring support k to 
that with k replaced by a rigid support. 
In the simple case of § 67 the foundation displacement 
was taken as proportional to the force because of 
its comparative rigidity. This cannot be done, however, 
where the foundation itself is relatively yielding and 
may introduce an additional resonant frequency. We 
have now to deal with a compound system. The floor of a 
building on which an elastically supported machine is 
mounted may constitute such a system, «Z 2 and repre- 
senting the equivalent mass and stiffness of the floor. In 
this case, e must be calculated in terms of foundation 
vibration amplitudes rather than transmitted forces, with 
and without the elastic suspension ki. 
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From the rules of § 31 

— miw")— 3’2 ^i = F 

y-ziki + ^2 ~ iiisor)— yd’i — 0 


(144) 


Sohung for j2, the foundation amplitude using the elastic 
suspension of stiffness 


3’2 = 


Fkt 


kih — (A’lWi -f- kini2 + k2mi)ccr + wi;?;2 


- (145) 

or 


When ki becomes relatively ver^'’ large, as with no 
elastic support between the vibrating body iui and floor 
mz, yz assumes the value 



F 

kz — {}Jii + ;»2) «- 


(146) 


yo 

By definition, the transmissibility = -=7^ == e 

3^ 2 


€ == 


kikz — k\(mi + Wz)o3~ 


kikz — {kiVii -}- kiviz + kziH\)oi' -b 7nin72W* 
Dividing the numerator and denominator each by kikz, 


(147) 


Vti + lllz 


03 - 
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JUlVlz 


03^ 


(14S) 


When the floor stiffness kz becomes very'^ large, this expres- 
sion reduces to 


1 mi , a>= _ (^2 

i 7- w- 1 - 

03c- 


(149) 


where w,. = y/kijiui. The transmissibility is thus seen 
to reduce to the form given by equation (140) derived pre- 
viously for the case of comparatively rigid foundation, 
as It should when kz is relatively large. 
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Expression (148) can be put in more convenient form 
if the two roots of the denominator are represented by 

o T , • 

oji* and and ; = ws'. In this case 

«/l + W2 


W3- 




6 ) 3 “ 


/ ur — W2“ — 

\ 031 “ / \ 032 “ / 


. . . (150) 


See reference 4. 

The constant 033 and one of the roots of the denomina- 
tor, 031 say, are both much larger than the other root 032 
in ordinary cases, and with a fairly stiff floor where 033 and 
031 are large compared with 03, the root 032 approaches 03 c = 
■\/kifmu so that expression (150) becomes 


1 




o o 

0)c“ — 03“* 


03c“ 


which is identical to (149) just preceding. Note that ai 
and 0)2 are the real resonant frequencies of the compound 
system (in radians per second), whereas 033 is not a real 
frequency but is the resonant frequency which the system 
would have if ki were a rigid connection. 

71. AppUcatiott to Single-phase Electric Motors. — In 
all single-phase motors the energ^"^ input must necessarily 
go through zero twice per cycle of the applied alternating 
electromotive force. This results in a double-frequency 
pulsating torque characteristic of all single-phase motors 
(reference 18), regardless of design except where special 
provision is made for energy storage as in the case of the 
capacitor type of motor. For small motors the armature 
current and consequently the periodic torque are rela- 
tively large at no load or small load. This torque rises and 
falls through an amplitude approximating the full-load 
torque in magnitude and at the magnetization frequency 
which is twice the alternating-current frequency. Thus 
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small, 60-cycle, single-phase motors when not elastically 
mounted correctly produce a 120-cycle hum which may be 
very objectionable. The hum from the average commer- 
cial washing machine may be transmitted all through the 
house where it is used. A correctly designed elastic 
mounting which will eliminate this hum can be applied 
at small cost. 

Figure 64 shows a mounting of theoretically correct 
design. Since the periodic torque responsible for the hum 
has a frequency of 120 cycles, tending to oscillate the stator 
about its axis at this fre- 
quency, the whole motor 
must be mounted so that 
it will oscillate about this 
axis at a natural frequency 
at least as low as -J of 
120 = 20 per second. 

Each of the end plates of 
the motor is connected to 
the base through a rubber 
ring fi.xed to an inner and 
an outer brass ring. The 
inner brass ring is at- 
tached to a boss on the 
end plate surrounding the 

shaft, the outer brass Fig. 64 . — Motor with Elastic Suspension, 
ring is fixed in a hole in 

a riser from the base plate. The result is that torque on 
the motor causes the rubber layer between the brass rings 
to shear a small amount. The shearing elasticity of this 
ring is such that the free torsional oscillation frequenc}'^ 
of the motor so mountedj is about 20 cycles. This type 
of mounting has the advantage of being very resistant to 
sideways displacement as produced by belt pull, and at the 
same time cuts out the 120 hum due to periodic torque. 

Other forms of elastic suspension can be successfully 
used for single-phase motors, but for best results the elas- 
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tic material should be placed as near to the shaft as is prac- 
ticable. If this is not done, and the mounting is made 
rigid enough to withstand belt pull, its torsional stiffness 
is very likely to be too great to be consistent with the 
torsional moment of inertia of the stator in producing 
a low natural frequenc3^ This greater stiffness is due 
to the longer moment arm when the support is at a greater 
distance from the stator axis. A support designed with- 
out calculating or measuring the particular natural fre- 
quency required will often afford little benefit. The 
majority of such elastic supports which the author has 
examined have been found to have too great torsional 
stiffness for good isolation. 

From § 40 it will be recalled that a torque acting on a 
mass free in space causes it to rotate about its center of 
gravity regardless of where the torque is applied. Thus 
a single-phase motor tends to oscillate about its center of 
gravity which usually is about on the stator axis. 'If the 
motor has other masses attached to it, its center of gravity 
may be considerably displaced. The elastic support 
should take account of the consequent displacement of its 
natural axis of oscillation in such cases, and be designed so 
that the axis of elastic torsional motion coincides with the 
natural axis of oscillation. If this is not done, transverse 
forces will be transmitted through the support which the 
support may be transversely too rigid to cushion. 

72. Application in Office Buildings and Hotels. — 
Modern office buildings and hotels must not be subjected 
to hums transmitted from electrical machinery. It is, there- 
fore, essential to design ventilating systems with great 
care to guard against noise from power-driven fans. For 
this reason it has been the custom to avoid the use of 
alternating-current motors in such applications because 
even the best alternating-current motors will hum some- 
what. Three-phase motors can be used successfully, how- 
ever, with suitable vibration-absorbing bases. If cork 
cushioning pads are used, 'as they often are, the motor 
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should be bolted to a hea\^’^ concrete base and the cork 
placed under this. If this is not done the cork pad may 
transmit considerable hum. To be effective, the cork sup- 
ports should be sufficiently thick and carrj' enough load to 
be well compressed. Without sufficient elastic compression 
(§9), the natural frequency of the support will not be 
lowered sufficiently to make the suspension effective. 

Suppose the motor has a 120-cycle hum, and also the 
upper harmonics of 240 and 360 cycles are present, which 
the best three-phase motors exhibit to some extent. An 
elastic suspension which gives an up-and-down frequency 
of 120 6 = 20 cycles should sufficiently screen the 

transmission of all these frequencies, provided that the 
five other natural frequencies of the suspension are so 
placed that they are not excited. 

From § 9 the compression of the elastic support neces- 
ssLTj to produce tliis frequency in an up-and-down direction 
can at once be calculated as follows: 

/ = 20 = = 0.0245 inch. (151) 

Thus a compression of about 25 mils will give the re- 
quired frequency, which from (140) gives a transmissi- 
bility of for 120 cycles and a much smaller transmissi- 
bility for the higher harmonics. Cork, however, is not 
a perfectly elastic material, and tests show that it should 
be compressed 3 or 4 times as much as an ideal spring to 
produce the required frequency, or as much as xV inch 
in the case just cited. (See reference 4.) If this pre- 
caution is taken, experience shows that it gives reasonably 
good results, though with the passage of time it tends to 
lose its elasticity somewhat under load. Considering its 
poor elastic qualities from the test standpoint it is sur- 
prising how well it works in some applications. 

From the standpoint of elastic suspension the purpose 
of the concrete mass is to insure sufficient compression of 
the cork base to give it a low enough up-and-down fre- 
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quency and at the same time to permit the use of enough 
cork to insure sufficient transverse rigidity to withstand 
belt pull. With a little practice a designer can easily 
determine the weight of the slab and thickness and area of 
the cork base necessary for the required degree of com- 
pression. 

High-grade rubber is elastically much superior to cork, 
nor does it lose its resilience appreciably with time, as is 
commonly believed, unless cheap grades are used. 

It is possible to design satisfactory rubber suspensions 
for fan motors and for other such applications which have 
the advantage of requiring no concrete mass. 

Expert knowledge is necessary, however; otherwise the 
support may be worse than useless. An important pre- 
caution in this type of suspension is to keep the natural 
frequencies of the support away from the operating speed 
of the motor, to avoid resonance in the operating range. 
If this range is considerable the design problem is corre- 
spondingly difficult. Every such suspension has six reso- 
ant jreqiiencies-. three possible angular vibrations, one 
about each of three perpendicular axes; and three linear 
vibrations, one in line with each of these axes. In other 
words, there are six degrees of freedom, each of which has 
its own resonant frequency (reference 17). 

For satisfactory results the elastic support should be 
designed by the motor manufacturer, with a knowledge of 
the operating range of the motor. 

It is difficult to design good steel spring suspensions of 
this type. Such springs do not offer the flexibility of 
design that rubber does in obtaining the right elasticity 
for every degree of freedom, and at the same time the 
necessary rigidity to carry belt pull. Furthermore, the 
small damping in rubber is a safeguard against vibration 
in passing through resonance. 

Figure 65 shows a motor in service mounted on a spe- 
cially designed rubber suspension. The rubber pads, 
placed in the metal boxes which protect them, are so 
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shaped that they have the correct stiffness in every direc- 
tion. For a further discussion of this subject sec reference 17, 



Fig. 65.- — Elastic Suspension Installation. 


73. Isolation from Vibration. — ^Thc reciprocal case from 
that of § 67 shown in Fig, 60 is illustrated in Fig. 66. 



Fig. 66. — Vibration. 
Transmission to Mass. 


Instead of finding how much motion or 
force is transmitted to the foundation 
with a given excitation of the mass ?;/, 
the question here is to find out how much 
motion the mass m receives through the 
suspension k when tlie foundation vi- 


brates through a fixed amplitude y 2 . 
The transmissibility is defined as the ratio of the ampli- 
tudes ji to y 2 , or € = yi/ 3 ’ 2 . From § 31 the vector equation 
for the forces on m is given by 


yi{k — muS^) — ky 2 = 0 
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Thus 


e = ^ = 




3’2 k — mor 


COc 


(152) 


which is exactly the same as equation (140) for the recip- 
rocal case. The expression for e when the spring contains 
damping is also the same for either case. 

Isolation of this character is important where a piece 
of apparatus must be screened from the effects of a vibrat- 
ing floor or building. Galvanometer suspensions and the 
balancing machine suspension taken up 
in § 63 are illustrations. 

74. Vehicle Suspension. — An im- 
portant practical application of this 
type of isolation is the spring suspen- 
sion of vehicles such as automobiles. 

Figure 67 shows a simplified diagram 
of the mechanical vibrating system of 
an automobile where the possibility of 
angular forward and backward pitch- 
ing is omitted to simplify the problem. The elastic constant 
^2 represents the resilience of the rubber tire, and ki the 
elastic constant of the springs of the vehicle. Setting 
down the vector equations according to the rules of § 31: 

yi{ki — yoki = 0 

y2{ki + k2 — 7??2W-) — yi^l — 3’3^2 = 0 





Fig 67. — Vibration 
Transmission to Com- 
pound System. 
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where and ar are the two roots of the denominator of 
equation (153) at which the system resonates and the 
transraissibility is indefinitely magnified. At points away 
from these resonances the transmissibility is small. 

Note that yi is the motion of the vehicle and ys repre- 
sents the amplitude imparted by bumps in the road. The 
higher frequency root, wi sa}^ corresponds to the resonance 
where mi and m 2 are vibrating towards and away from 
each other out of phase. “Washboard” roads causing 
the wheels to resonate up and down against the stops 
on the A^ehicle body illustrate this resonance. The lower 
resonance takes place when a car goes over much longer 
waves in the road, making the whole car body resonate up 
and do^ra. 

The damping produced by shock absorbers, not taken 
into account in the analysis presented, is essential in reduc- 
ing the amplitudes of these two resonant frequencies, since 
every car must sometimes be subjected to one or the other 
of these two resonant frequencies on an uneven road. 
The function of shock absorbers is thus to damp resonant 
oscillations. So far as reducing shocks between resonant 
frequencies is concerned, they actually increase the trans- 
missibility slightly, as explained for a simple case in § 69. 

In general, the softer the springs the better, so far as is 
consistent with the stability of the car. This subject 
offers a field for further study. There is little doubt that, 
with careful study from the standpoint of these funda- 
mentals, more can be accomplished in this direction than 
is realized by manufacturers. 
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GENERAL THEORY OF VIBRATION DAMPING 

75. Vibration Damping and Logarithmic Decrement. — 
The study of vibration damping is a whole field in itself 
which is of decided interest and Importance in vibration 
work. 

If a system of one degree of freedom with damping 
proportional to velocity, like that of § 28, vibrates freely 
with no stimulating force F, the sustained vibration term 
drops out and there remains only the expression 

cos cooi (155) 

Figure 68 is a 
graphical representa- 
tion of this equation. 
The best way to specify 

^ the damping in a vi- 

^ brating system is by 
means of its loga- 
ritlimic decrement b, 
Fig. 68.— Decaying Vibration. often spoken of simply 

as decrement. 

The logarithmic decrement is defined as the natural 
logarithm of the ratio of two successive vibration amplitudes 
during the free vibration. At P in the vibration decre- 
ment curve shown, 5 = log {y\/y^. 

From (155), when / = 0, yi = yo. 

One cycle later, t — T, the period, and y^. — yoe~“^. 

5 = log^ = loge“^ = aT . . (156) 

^'2 Jo 

where /o = natural frequency = 1/T. 
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In this case, since both a and /o are constants which are 
independent of amplitude, 5 = a constant at all times 
during the \dbration decay. Thus a free vibration whose 
amplitude falls off exponentially with time like that given 
by equation (155) has a constant logarithmic decrement. 

In practice, 5 is seldom a constant, but either falls off 
or increases as the amplitude of the free vibration decreases. 
A velocity damping, however, where the frictional force 
is proportional to velocity (friction force = hdy/dt) always 
yields an exponential amplitude decrement cur^'e like that 
of equation (155) with a constant logarithmic decrement. 

Decaying electrical oscillations also fall off exponen- 
tially since electrical resistance acts as a velocity damping. 
The damping force in such an electrical circuit is a counter 
voltage of magnitude R{dg/dt), where R = electrical resist- 
ance analogous to the friction constant 
b, and q is the charge displacement 
analogous to the mechanical displace- 
ment y. 

Although a velocity damping pro- 
duces a logarithmic amplitude decre- 
ment of a free vibration this type of o9.-Sin,plc DamM 
friction is not the only kind which will Vibration, 

do this. The condition which has to 
be fulfilled for logarithmic amplitude decrement is simply 
that the dissipation of energy' per cycle be proportional 
to amplitude squared, which condition velocity damping 
satisfies. This condition is also satisfied by a damping 
proportional to amplitude squared, which is independent 
of frequency. 

76. Dissipation per Cycle and Loop Area. — The pres- 
ence of damping in the elastic support of a weight IF as 
shown in Fig. 69 causes the load deflection cuiwm to 
depart from the ideal straight line AB, Fig. 70, and assume 

a curve which describes the loop ACBD for a complete 
cycle. 

Just as in the case of the magnetic hysteresis loop, the 
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area of this loop is a measure of the energy dissipated during 
the cycle of motion. This is easily shown to be true from 

the principle of work. The 
loop is described in the 
direction shown by the 
arrows. The force is measured 
upwards from the horizontal 
line of zero force at the 
bottom of the diagram and 
rises and falls above the 
mean value fo during the 
cycle of motion. The motion 
ismeasured horizontally from 
the vertical coordinate axis 
which corresponds to the rest 
position, the downward ex- 
tension of the spring corre- 
sponding to increasing force 
being taken positive, and up- 
being taken as negative displacement. 



70. — Stress - Strain 
Damped Vibration. 


ward relaxation 
Now 


Work 


=fjdy 


(Note that / = force and not frequency as used here.) 
During the motion ACB this work input is positive since 
the work is done upon the spring and is equal to the area 
EACBF; during the motion BCA the work is supplied by 
the spring and is negative and equal to the area EADBF. 
Thus during each cycle the work input exceeds the output 
by the loop area. This loop area is thus equal to the dis- 
sipation per cycle in terms of the force and displacement 
units used. If / is in pounds and y in inches, the loop area 
measures the dissipation during the cycle in inch-pounds. 
The dissipation during one cycle will now be calculated. 
From § 31 the vector equation of forces is 

yi{k — mur +jbu))= F ... . (157) 
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jf y = sin a'/, all terms in the left-hand member are 
sine functions except the j term which is 90° ahead of the 
sine terms and is thus a cosine term. Thus, dropping the 
vector notation 

— war) sin a?/ -}- 3’i^w cos ut = f . . (158) 


where f is the instantaneous sinusoidal force, whose phase 
is determined by the left-hand member of the equation. 

Work per cycle = fdy over a complete cycle 

jsdy = if [yi(/.’ — mur) sin ut 3 'i boi cos w/jri cos ut d {ut) 
since 

dy = 3 'i cos uiid (uf) 


Integrating over a cycle, the sin cos term is zero. Therefore 

^ fdy = 2 yrbu cos-o.'l d{o}l) = -yrbeo (159) 

This is tlie work done per c\'clc by the force / and is there- 
fore the dissipation per cycle due to the friction. 

In this analysis, / is measured from the center of the 
loop instead of from the base line, thus leaving out the con- 
stant gravitational force component. The gravitational 
force /o jdelds no work in a complete loop which ends at 
the same point it started from since work done b}'^ graA'ita- 
tion due to an equal up-and-down motion is equal and 
opposite in sign, when the mass is brought back to its orig- 
inal position. 

The expression Tryr&co equals the loop area in units of 
/ Md y. The loop can be shown to be a perfect ellipse for 
this type of damping. Its ciirs'e is given by the equations 

f = 3'i(^ — mor) sin cot -|- yibco cos co/1 
y = 3 ’i sin co/ J ^ 
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Eliminating the parameter wt by substituting, squaring, 
and using the relation sin- at -f- cos- w/ = 1, there is 
obtained after simplifying, 

[(ife - + b‘^a-]y^ - 2{k - wco2) yf + p = (161) 

This is an equation of the second degree in y and / and 
must thus be a conic section, and since it is a closed cur\'^e 
it is an ellipse. The yf term is present because the ellipse 
is tilted with respect to the coordinate a.xes. At reso- 
nance, k — = 0, and equation (161) becomes 



P 

yiW-oP 


= 1 . 


(162) 


When the elastic force ky and mass reaction —moi^y balance, 
as at resonance, the only remaining force is that due to 
damping, and the ellipse (Fig. 69) becomes horizontal, 
as seen from equation (162), with a semi major axis equal 
to the amplitude of deflection y\ and the semi minor a.\is 
equal to the frictional force amplitude yj6w. The area 
of this ellipse in terms of its semi axes equals iryi^ba, which 
checks equation (159). 

, Note that below resonance the ellipse twists counter- 
''lockwise from the horizontal position, and above reso- 
nance it twists clockwise. It is a twist in which every 
vertical element retains the same breadth but whose 
heighv depends upon w. Thus the total area is propor- 
tional 1 1 frequency, as it must be with velocity damping. 

Whv y and / are replaced by strain e and stress cr at 
some pa T in an elastic medium supposed to contain a 
velocity 't iping force T}{de/di) corresponding to bidyjdt) 
of the pT^ us problem, the loop area is given by an 
expressions t (159) 

\ Area = irem-rja (163) 

where £„ represents the amplitude of the sinusoidal 
strain e, and rj the coefficient of normal viscosity. In this 
case the dissipation is'for unit volume of the medium per 
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strain cycle at the point under consideration. For another 
deri^^tion of the expression see § 79. 

77. Energy Expression for Decrement. — ^An important 
formula null now be derived which connects the logarithmic 
decrement Snath the dissipation per vibration cycle just dis- 
cussed, which will be called ATF, meaning a small amount 
of work. 

Let TT^ equal the total vibrational energy’’ of some vibrat- 
ing system vibrating in a given mode. It will be recalled 
from § 12 that, when a system vibrates in a given mode, 
the Aubration of all of its particles is isochronous, that is, 
they \'ibrate all togetlier in the same period, all reaching 
their maximum and minimum displacements at the same 
instant. Then IF = Ky-, where K is a. constant depend- 
ing on the system, the mode of vibration, and where y 
is measured. To be specific, take the cantilever shown in 



Fig. 71. — Vibrating Cantilever. 


Fig. 7 1 vibrating in its gravest mode, and let y be measured 
at P. It is evident tliat K is different from the elastic 
constant k, and is a general proportionality constant. The 
e.\pression W = Ky- says that the vibrational energ^^ is 
proportional to the amplitude squared, the magnitude 
of K being such as to give the correct IF for the case 
considered ^ith y measured at P. This expression follows 

from Hooke’s Law tliat stress is directly proportional to 
strain. 


^ During a decaying vibration the amplitude at P grad- 
ually falls off. Assume that during one cycle it falls off 
from y to y — Ay. The dissipation per cycle ATF must 
be the cause of the amplitude decrease Ay. 

Since 

IF = A"y= 

IF - ATF = K{y - Ay)2 

= K{y- — 2yAy -f- Ay-) 
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therefore AW = K-2yAy, neglecting Aj?-, which is small 
compared with 2yAy. 

Therefore 


ATT K-2yAy Ay 
2TT “ 2Kf ~ T 


(164) 


Bj'’ definition (§ 75), 

; , tog y+ ^ y , ^ 


(165) 


where the small terms in higher powers of Ay/y are dropped 
because theji- are relatively small. This result is obtained 
by expanding the logarithm into its series. 

Comparing (164) and (165) it is seen that 


ATT _ 
2T4^ ^ 


(166) 


Thus 5 can be expressed in three ways. 


( 1 ) 

( 2 ) 

(3) 


5=j^(§ 75) 

y 

2 IT 


(167) 


The first is an exact expression. 

The second and third are approximate, but the approxi- 
mation is good for practical work where Ay is small com- 
pared with y and thus A W is small compared with IT. 

The second expression Ay/y affords a quick method 
of estimating 6 from the graph of a decaying vibration 
like that of Fig. 68. 

The third expression is extremely useful in practical 
work in vibration damping as giving at once the energy 
dissipation produced when 5 is known, or vice versa. If 
ATT can be estimated for a damping device, 5 can be imme- 
diately found, giving a quantitative estimate of the effec- 
tiveness of the damping means. 



ENERGY FORMULA 


117 


The great advantage of- expressing the damping of a 
vibrating system in terms of 5 is that 5 is a dimensionless 
quantity, being a numerical ratio, and is thus the same 
whatever system of units be used to express energy and 
length. It is easily measured. It is readily found from 
the results of other investigators of vibration damping 
phenomena, however the results be set forth. By means 
of the energy expression when 5 is known, AW is easily 
found, and vice versa. 

The formulae (167) have great generality in that they 
apply at any instant, that is, at any point on an amplitude 
time decrement curve like that of Fig. 68. They apply to 
any system vibrating in any mode whose vibration decay 
can be represented by such a decrement curv'^e. Therefore 
they apply regardless of how damping be produced, Avhether 
by viscous friction, externally or within the spring material, 
whether by coulomb friction or solid friction, and of course 
regardless of how the damping depends upon tlie amplitude 
of motion, or what be the shape of the “hysteresis” loop 
area during the vibration cycle. See reference 19. 

78. Illustrations. — ^The application of the formulae for 
decrement just derived will be illustrated by a few simple 
cases. 

(1) Comparing the first of equations (167) with the 
transient term of the solution of equation (44) of § 28, 


where m = W/g. 



(168) 



Therefore 


b 

2mfo 


But/o = 




1 IK 

2tv yin 


nearly 


5 


irb 

y/ km 


( 169 ) 


which is a constant for this system. 
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Note that the vibrating sj'stem of § 27 for which tliis 
value of 5 was found is the same as that analyzed in § 75 
and shown in Fig. 69. For this system it Avas found that 
AW = iry'^bw (dropping the subscript of 3 '). Also for this 
system W = ^ky~. Note that in this e.\pression K of 
§ 77 = 

From the previous section 5 = ATF/2TF. 

Substituting the values of W and ATF, 


5 = 


Trby~o3 

ky^~ 


irbw 

k 




m 


since, as before, <0 = Wc = y/kini very nearly, because 5 
is expressed in terms of w for a free vibration. Note that 

(170) checks with (169) derived just 
/A ^ I 1 previously. 

(2) Suppose a mass m to be vi- 

^ brating back and forth and to be 

Coulomb Friction Damping, gradually brought to rest by sliding 

or “Coulomb” friction as indicated 



in Fig. 72. 
plete cycle, 

also 

Therefore 


If the frictional force is F, considering a com- 
ATF = 4Fy 
W = Py2 


S = 


AW 

2W 


4Fy 

ky"^ 


ky 


(171) 


Note that 5 is not a constant but increases as y decreases. 

In this case the amplitude falls off according to a 
straight line as shown in Fig. 73. This can be shoAvn as 
follows. 

For any vibration, 

W = Ky^ (§ 77) 
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But 


dW 


= — LW X frequency = — AWf 
dt 

Substituting and solving for dy/dt, 
dy — AWf 


dt 


2Ky 


(172) 



Fig. 73. — Decaying Vibration with 
Coulomb Friction Damping. 


Thus if AW is linear in y, 
dy/dt = constant; that is, the 
amplitude falls off with time 
according to a straight line. 

Equation (172) gives an ex- 
ponential decrease of ampli- 
tude if AW depends on y^ as 
previously noted, and if AW 
depends on y^ the ampli- 
tude time decrement cur\’^e is hyperbolic, equation (172) 
giving the slope of the decrement curve in any case where 
ate can be expressed as a function of y. 

79. Viscous Damping in Spring Material. — Suppose a 
mass m to be fixed to a reed which is 
rigidly attached at the lower end as 
shown in Fig. 74. Assume that all the 
friction present in the system causing 
vibration decay resides in the material of 
the reed itself and that this friction is 
viscous in character. Every elementary 
volume of the reed material goes through 
cyclical strains which for one vibration 
cycle are given by e = sin iot, neglecting 

Fig. 74.— Vibrating falling off of amplitude in one 

Mass on Reed. cycle. 

If the internal friction is of viscous 
character a resisting force is present which may be specified 
by the expression / = rj^de/dt). This expression says that 
if unit volume of the material of the reed is strained with 
unit velocity, the frictional resistance/ = ij. The constant 


m. 
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•jj is the viscosity coefficient as defined by / = riidt/dt) and 
is perfectly analogous to b used in § 76 for the system of 
one degree of freedom. In this case a point in the solid 
material of the reed is considered, and is expanded to occupy 
unit volume by the imagination; in the case cited where b 
is used the whole system is considered as a unit. 

As explained in § 76, if AH-^i = dissipation per tinit 
volume 

APFi = TreJ^rjo) .... (173) 

ATI^i is given the subscript because it represents dissipa- 
tion per cycle per unit volume instead of for the whole 
vibrating system. 

For the same unit volume the maximum strain energy 
stored in it per unit volume 

IFl 

Since (t„ = Etm (§ 14), 

Wi = iEeJ 

where Cm and e„ are maximum values of stress and strain 
respectively at that point. 

Therefore for unit volume 

^ ATFi ir-qeja Trijic 


If this single unit cube were to vibrate by itself, (174) 
gives its decrement for viscous friction. To find 5 for the 
whole reed, however, the value of AW and of W for the 
whole volume V of the reed must be used. For the entire 
reed material. 


ATT _ /// Tr7jem“a«/F 

“ JJ/ EiJdV 

TTVU) III TTVU) 


( 175 ) 
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which is the same as (174) for the single unit volume. It 
can be seen that this must be so because, for every unit 
volume, the ratio LW jlW is a constant for any type of 
vibration, for both depend upon The integrals in 

(175) cancel off since they are identical. 

The total vibrational energy^ lies in the reed as poten- 
tial energy of bending at maximum deflection. All the 
dissipation also takes place in the reed. Every unit volume 
of the reed contributes dissipation and vibrational energy 
to the system in the same ratio, regardless of strain ampli- 
tude. 

Note that, for the whole system vibrating, 8 is much 
less than it would be for unit volume vibrating alone, for 
although 8 is given by the same e.xpression, ojo for unit 
volume is much higher than is coo for the system as a whole. 
Although this is true of viscous damping, there are other 
types of damping where 5 does not depend on co and so 
is the same regardless of frequency. 

It should be pointed out that 8 is by definition based 
upon the natural frequency of the vibrating system con- 
sidered and so is a physical constant of that system. Thus 
in the final result w is always equal to coo, corresponding to 
the natural frequency /o. In some of the previous cases 
coc has been used because of its close appro.ximation to 
(OQ (see § 27). 

80. Solid Damping in an Elastic Material. — It has 
been shown beyond doubt that 6 for practically all solid 
materials shows no such dependence upon frequency as 
given by (175). Common observation points toward this 
conclusion, since if this were not so the higher harmonics 
of a vibrating body, such as a bar, for example, would be 
immediately damped out because of their higher 8, leaving 
only the lower harmonics. A bell with a tone of compara- 
tively high frequency would not ring at all, but act almost 
as dead as though made of putty. The viscous law of 
friction does not apply to metals nor to other solid materials 
as has been assumed by many investigators. 
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Experiments by the author * and by others show that 

AT'Fi = . . . . (176) 

approximately when the stress amplitude is considerably 
below the elastic limit. This expression is accurate so 
far as independence of frequency is concerned. The 
dependence on e„“ must be considered as a rather rough 
approximation. If tr„ approaches the elastic limit, €„ 
usually has an e.xponent greater than 2. For practical 
work, however, where the vibrational stresses are mostly 
far below the elastic limit, expression (176) will suffice. 
This expression has an advantage from the standpoint of 
analysis that great simplification results compared with a 
law where is raised to a power other than 2. 

Note particularly that expression (176) is independent 
of frequency. The constant ^ is merely a proportionality 
constant or solid friction constant, and depends upon the 
system of units employed as does E. 

This type of friction has been called solid friction as 
being the type of internal friction characteristic of most 
solid materials. 

For this law of dissipation. 


EWi = 

and as in the previous section, for unit volume, 

Wi = lEeJ 


therefore 


S = 


ATFi 


2TFi 




. (177) 


It also can be shown as in the previous section that for the 
whole elastic mass 

d = ^E 

In this case, however, every harmonic in the vibrating 
system has the same 5, regardless of its frequency. The 
higher harmonics die out most rapidly because they per- 


* See reference 20. 
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form a larger number of vibrations in a given time than 
the lower ones although 5 is the same for all harmonics. 
Here also, we know from common observation that the 
lower harmonics of a vibrating bar persist for a longer 
time than the higher harmonics, although the latter do 
show considerable persistence. 

This is an important conclusion, since when 5 is once 
measured for a given elastic material it is known for any 
vibrating system whose elastic energy and damping lie in 
that material, regardless of the frequency of the system. 

The potential energy of a vibrating system is frequently 
partially or wholly derived from gravitation, a gravity 
pendulum being an example of the latter. In this case 
the damping friction which brings it to rest comes from 
pivot friction and air friction. 
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SHAFT WHIRLING DUE TO INTERNAL FRICTION AND 
TO OIL ACTION IN JOURNAL BEARINGS 

81. Intemal Friction and Shaft Whirling. — It seems 
rather surprising that intemal friction within the material 
or structure of a revolving shaft or rotor should produce 
whirling. Ordinarily friction of this character in an 
elastic body causes vibrations to fall off and disappear. 
It is nevertheless true that friction in a rotor may cause 
it to build up a whirl under the right conditions. 

Figure 75 (a) shows a vertical shaft carrying a fly- 
wheel deflected to the right, and (b) shows a section of 
the shaft at CD near its middle point. The shaft and 
flywheel are supposed to be whirling about the center line 
of the bearings AB, the shaft center at the section CD 
describing the circular path in the direction of the arrow,, 
shown dotted in Fig. 75 (p). Since the elastic stiffness of 
the shaft balances the centrifugal force of the flywheel 
during the whirling motion, this whirl must take place at 
the angular velocity coc = V k/m, where k is the elastic 
constant of the shaft at the point where the flywheel is 
attached, and vi is the mass of the latter. The mass of 
the shaft is supposed to be small compared with that of the 
fljnvheel which it carries. Note that the shaft may revolve 
about its geometrical center 5 in either direction at any 
speed while the whirl takes place of fixed angular velocity 
lOc. The latter motion may be compared to that of a 
planet about the sun, and the former to the revolution of 
the planet about its own axis. The two motions are inde- 
pendent of each other. 

The whirl Wc is resonant in character and is easily built 
up. If the whirl took place at any other speed it would 
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have to be forced and thus be exactly in step with some 
periodic exciting force. In this case no force is supposed 
to exist other than the steady disturbing force Ff pro- 
duced by internal friction, which is always approximately 
tangent to the path of whirl. Since Ff acts in line with 
the direction of whirl, it does a definite amount of work 
per cycle in building up and sustaining the whirling motion. 
\\dien the rate of energ>^ dissipation just balances the 
rate of supply, the whirl will build up no further. 


C- 



Fig. 75. — Shaft Whirl Excited by Internal Friction. 


The force Ff is produced by internal friction within the 
rotor, and is in the positive direction of whirl only when the 
rotational speed of the shaft, w, exceeds its whirling speed 
Wc If w is less than coc the force Ff reverses and whirling 
is damped out. 

Consider a section of the shaft as shown in Fig. 75 {h). 
The right-hand half of the section of the shaft is under 
elastic tension Te. The left-hand half is under elastic 
compression Ce. As a fiber, whose section is shown at P, 
say, travels over the upper half of the section towards the 
region of elastic compression it is being shortened at a 
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certain rate. The internal friction present exerts a force 
of reaction to this shortening, so that a component of 
compressional force of small magnitude is exerted on this 
fiber. In like manner, a component of tension is required 
to overcome friction during elongation as the shaft fibers 
traverse the bottom half of the section. The resultant 
effect of this action is a small component of compression 
Cf in the upper half of the shaft section, and one of ten- 
sion Tr in the lower half of the section, due to internal 
friction. The elastic tension and compression Te and Ce 
due to bending produces a strong reaction Fe towards 0, 
balancing the centrifugal force of the flyAvheel m due to 
the whirl Wc The frictional components Tp and Cf in 
like manner produce a much smaller reaction component 
Ff along the path of whirl, which can be balanced only by 
forces resisting the whirl. If these forces of resistance 
are relatively small, the whirl will build up. 

This reasoning is general, and will apply to any kind 
of damping friction whatsoever in a rotor, whether it be in 
the shaft material or due to friction between parts of the 
rotor as its sides are alternately extended and compressed 
during its whirling motion as just described. 

If the angular velocity of the shaft w is the same as 
the whirling speed no friction will result as all the ele- 
ments of the rotor remain strained the same amount with- 
out the cyclic change of strain depending on the difference 
between w and Wc The shaft then would whirl with the 
same side out all the time, as though it had a permanent 
set. Thus for Ff to be positive the angular velocity of the 
shaft w must exceed Wc, the velocity of whirl. See refer- 
ence 21. 

82. Occurrence in Practice. — The chief danger from 
this type of whirling is that because of its rarity it is not 
recognized and understood in the few cases where it does 
occur. It is easy to produce such a whirl in a model 
shaft with a flyAvheel if one or two rings of axial length as 
great as the diameter of the shaft are shrunk on near the 
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middle of the span. Friction arises during deflection from 
a slight axial slipping of the shaft surface inside of the rings 
which is sufficient to cause a vigorous whirl when the shaft 
is run up above its critical speed. Tests show that the 
internal friction in the steel of a shaft is not sufficient to 
build up such a whirl. Built-up rotors which operate 
above the critical speed, however, must be watched, 
and should be so constructed to avoid as far as possible 
friction between parts when de- 
flected. The impellers of com- 
pressors and disks of turbines 
should have narrow axial contacts 
vdth the shaft where they are 
shrunk on. Long, flexible lami- 
nated rotors which vibrate badly 
regardless of balance are open to 
suspicion. At least one case of a 
bad shaft failure probably due to 
this cause has come to the author’s 
attention. 

83. The Oil Whip. — It has 

been found from practice that the 

action of oil in journal bearings is 

a more common cause of shaft 

whirling than internal friction in 

the rotor. This type of whirling 

will be taken up at this point be- 76.-Shaft Whirl Excited 
, ^ ^ . by Oil Action, 

cause ol its similarity to internal 

friction whirling, although it is not a damping or friction 
phenomenon. 

If a small model shaft (Fig. 76) about 1 inch in diame- 



ter carrying a fljnvheel is placed on end with a ball bearing 
B at the bottom and a journal bearing A at the top with a 
small clearance which is well supplied with oil, it will be 
found, on slowly revolving the shaft, that the axis of the 
journal does not remain stationary, but moves in a small 
circle. This can be observed very well by machining the 
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point P on the upper end of the shaft, and noting its 
motion under a microscope while illuminated. The jour- 
nal describes a circular motion in the bearing of amplitude 
depending upon the clearance. Furthermore, it will be 
noted that the angular velocity of this circular motion of 
the journal axis is exactly one-half of that of the shaft and 
flywheel. This motion is evidently produced by the action 
of the oil film, since the motion ceases when the oil supply 
is cut off. It can be explained by the action of the wedge- 
shaped oil film produced by the eccentric position of the 
journal. The oil film acts as a wedge and forces the 
journal around in the bearing in a forward direction at the 
speed with which the wedge of oil moves forward. This 
forward speed equals the average speed with which the oil 
wedge is dragged forward by the journal surface. One 
surface of the oil film is held stationary by the bearing 
while the other is dragged along at the journal speed, the 
average speed thus being about half the journal speed, 
just as the forward speed of the balls in a ball bearing is 
about half the journal speed. 

For slow rotation the wedging action of the oil is not 
strong, and the journal performs the circular motion 
only when almost balanced on its end in the vertical 
position. 

If the speed of the shaft is raised up to slightly above 
twice its critical speed it will then begin to whip, provided 
the bearings are both solidly supported so as to avoid the 
dissipation of vibrational energy. This building up of a 
whip is evidently due to the angular velocity of the wedge 
of oil falling in step with the critical speed of the shaft, so 
that a resonant vibration is established. When the shaft 
once starts to resonate at its critical speed, the oil wedge 
is held down to this speed even if the rate of revolution 
of the shaft is much higher than twice its critical speed. 
The driving force back of the wedge of oil increases, but 
this merely serves to stimulate the resonant whip of the 
shaft more strongly, while the velocity of the oil wedge is 
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held down to the critical speed frequency at which whip- 
ping takes place. 

The tendency for this type of whipping is quite strong 
at rotational speeds of 25 per cent or more above twice the 
critical, so that a rotor in the horizontal position may be 
subject to this type of whipping, in spite of the force of 
gravitation which tends to hold the journal down on one 
side of the bearing. As soon as an oscillation is started at 
the resonant frequency of the rotor, the building-up action 
of the oil wedge becomes active. 

That oil action is the cause of such whipping is shown 
by either cutting off the oil supply or reducing the speed 
of the rotor below twice the critical, since in either case the 
whipping ceases at once. The phenomenon is easily shown 
on a small model. 

It has been found that increasing the bearing pressure 
by decreasing its area decreases the whipping tendency 
from oil action, but this is not a certain remedy. It is 
interesting that twisting the bearing a.\is out of line from 
the journal axis stops the whipping almost invariably. 
This method is uncertain in practice, as such lack of align- 
ment is difficult to maintain and is apt to be soon lost 
during operation. 

Another way to stop such whipping is through the use 
of spring supports with friction damping, but such sup- 
ports as have been found effective are expensive and 
require considerable attention. The best cure is to avoid 
running journal sets over twice their first critical speed. 
See reference 22. 
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MEASUREMENT OF DAMPING CONSTANTS OF SOLID 

MATERIALS 

84. Revolving Rod Method of Measuring Internal 
Friction. — One reason for discussing in § 81 the theory of 
shaft whipping as produced by internal friction is that the 
measurement of the disturbing force Ff gives a quick and 
direct method of measuring 6 for the material of which 
the shaft is made. 



Fig. 77. — Revolving Rod Test for Internal Friction. 


Instead of a vertical shaft, consider one which is hori- 
zontal, and for convenience, one which is overhanging as 
shown in Fig. 77. 

If the overhanging shaft or rod CB is supported in the 
bearings C and D and driven through the pulley P, and 
at the same time its end is deflected downwards by the 
force F, produced by the weight and pan carried on another 
bearing at B, it will be found that, on revolving this 
deflected rod, it will also deflect a small amount sideways 
to the right or left, depending upon the direction of rota- 
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tion. Figure 78 shows a section at AB. The end of the 
shaft is deflected downwards by the amount h, and swings 
to the right through the angle <j). This angular deflection 
is caused by the internal friction in the shaft, and the 
magnitude of Ff (§ 81), which causes this deflection, 
is thus equal to Fsin^, while the elastic component 
Fe = F cos cf). 

Let us apply the formula B = AW/2W to this system. 

The torque required to produce the deflection (f> against 


the gravitational force / = Fph = Fh 
sin 4>. Thus the total work done on 
the rod per cycle in overcoming its in- 
ternal friction equals 
AWr = IttIiF sin cf, ... . (178) 

where AWn is the total internal friction 
Avork per cycle of revolution. Now if 
the gravitational force Avere eliminated 
and the mass m Avere vibrated up and 
doAvn through the same amplitude h, it 
can be shoA\m by an integration that 


ATT 


ATT;; 

2 


(179) 


In Avords, expression (179) says that 
half as much energy is dissipated per 
cycle in a flexural vibration as is dissi- 



pated per cycle in revolving the rod AA'hile 
deflected, the amplitude h being the same 
in each case. 


Fig. 7S. — Transverse 
Displacement of End 
of Rod. 


This foIIoAvs from the laAv of solid friction giA^en hy 
(176), § 80, AA^here dissipation depends upon stress squared, 
and may be shown as folloAvs. 

For the revolving rod, see Fig. 79, the fiber at P, of 
cross-sectional area rdddr, goes through a stress cycle pro- 
portional to r Avhile for the same rod Avith in vibrating up 
and doAvn the stress cycle of this fiber is proportional to 
r sin 0, Avhere AB is the neutral axis of the cross-section 
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for the vibration. Thus on the basis of § SO the dissipa- 
tion per cycle per unit length of fiber at P is kr- X rdddr 
for the revolving rod and kr- sin- 6 X rdddr for the vibrat- 
ing rod, where k is a proportionality constant, the same in 

each case, since h, Fig. 78, is taken 
the same for rotation as for vibration. 

For unit length of rod, these ex- 
pressions are integrated over the area 
of cross-section and will be found to 
differ by the factor 2. Assuming 

ra%ibra7ng''anT'RSting ^hat the deflection curvc of the rod 
Rod. is the same when it vibrates as when 

deflected downwards during rotation, 
the dissipation per cycle for the entire rod differs by the 
factor 2 for the two cases so that 



AWr = 2AW 

Therefore from equation (179) 

AW = irliF sin 


The potential energy W due to the elastic force com' 
ponent Fcos<l) (Fig. 78) = ^JiFcostf). Therefore 


d 


AW 

2W 


irhF sin <f> 
hF cos (j) 


= TT tan (j> 


(180) 


85. Damping Constants for Various Solids. — ^Thus by 
carefully measuring <{>, 5 can at once be found from equa- 
tion (180). Its value is approximate except for the case 
where AW depends strictly on stress amplitude squared, 
in which case AW = AWr {2 exactly. 

It is easily established that 6 and thus AW does not 
depend upon frequency since is entirely independent of 
the speed of rotation of the rod over a wide range, regard- 
less of the load 7n on its end. By this method many 
metals and other materials have been found to exhibit this 
characteristic, so that it may be said that an outstanding 
characteristic of solid friction is that the dissipation per 
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vibration cycle is independent of frequency as given by 
(176) in § 80, all the way from less than one to many 

■^ousands of cycles per second. 

By varying the force F the angle (f> in general shows 
some variation, usually increasing with load. For a con- 
stant 5, (j) must be constant also as is seen from expression 
(180). 

In the following table is given a series of values of 5 for 
different materials calculated by this method. These 
values are for the most part approximate as they depend on 
an average value of tan 6 for the material in question 
except in special cases where tan (f> is practically a constant 
for all loads and AW thus depends exactly on amplitude 
squared. Regardless of this approximation, these values 
of 5 approximate the truth sufficiently to be very useful in 
practical work. 


Table II 


Material 

Tan <}> 

S 

Rubber 90% pure * 


0.26 

Celluloid 

0144 

.045 

Tin, swaged 

0407 

.129 

Maple wood 

00669 

.022 

Zinc, swaged 

00630 

020 

Glass 

00205 

0064 

Aluminum, cold-rolled 

.00107 

.0034 

Brass, cold-rolled 

.00154 

.0048 

Copper, cold-rolled 

.00159 

.0050 

Tung:sten, s\vaged 

.00524 

.0165 

Swedish iron, annealed 

.00250 

.0079 

Phosphor bronze, annealed 

.00101 

.0032 

Mild steel, cold-rolled 

.00157 

.0049 

Molybdenum, swaged 

.00219 

.0069 

Nickel, cold-rolled 

.00102 

.0032 

Nickel steel, 31% siraged 

.00073 

.0023 

Monel, cold-rolled 

.000454 

.00143 

Phosphor bronze, cold-rolled to maximum hardness. 

.000117 

.00037 


* Vibration test. 
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VIBRATION DAMPING IN TURBINE BUCKETS 


86. Vibration Amplitude from Energy Input. — Suppose 
that a system of one degree of freedom with the elastic con- 
stant k and mass m which is known to have a decrement 5 is 
acted upon by a known pulsating magnetic force. It is 
required to find out the amplitude of vibration at resonance. 
A pulsating magnetic force can be very approximately 
represented by the expression 

/ = I? (1 -f cos w/) .... (181) 

where Jo is the total amplitude of the pulsating magnetic 
pull. Let the vibratory motion produced by this force be 
given by 

y = yo sin (wf -p a) . . . . (182) 


where a is the phase angle. 

The work done per cycle equals 

m =^Jdy=2jj^ {\+coswt)cos{oit + a)d{a) (183) 


At resonance, dy is always positive while / is positive, and 
vice versa, so that a = 0. Thus integrating (183) for a = 0, 


But 


ATV = 

ATT _ ir/oyo _ tt/q 
2 IT Ikyd^ 2kyQ 


W = 

134 


since 


. ( 184 ) 



VIBRATION OF TURBINE BUCKETS 


135 ' 


Solving for yo, 



(185) 


which is the -required relation. 

This expression gives the maximum possible amplitude 
obtainable in a vibrating systeni with decrement 5 from 
the pulsating force (184) acting upon the 
system at its resonant frequency. 

87. Vibration Amplitude of Turbine 
Buckets. — It is sometimes required to find 
the approximate vibration amplitude of a 
reed-like body, such as a turbine bucket, 
which arises from a periodic force of the 
type given in the previous section, the 
damping constant S being known approxi- 
mately. The maximum periodic stress at 
the root of the bucket can then be calcu- 
lated. 

Consider a bucket of length 1. 

At resonance the energy input per cycle 
equals the dissipation per cycle = ATT. 

If X is measured along the bucket and 
u IS the vibration amplitude at any point Vibration 

and/o is the maximum value per unit length from PeriodicForce. 
of the pulsating force, / = /o/2 (1 -f cos ut) 
distributed uniformly along the bucket, at some point P, 
the input per unit length is 

= (§ 86 ) 



Therefore the total input per cycle 

ATT = ^ fudx . . . 

Z Jq 

The dissipation also = ATT = 2TT'^5. 

From Fig. 80 



• . ( 186 ) 



136 


DAMPING IN TURBINE BUCKETS 


where w = weight per unit length. Thus dissipation equals 

Word 


2W8 = 


S 


- f u'^dx 

-Jq 


• (1S7) 


Equating input to dissipation and rearranging, 

f\fldx 


7r/og 

2w03~8 


f\idx 

•/o 


= a«i 


. (188) 


where a = ratio of the integrals for the case where U\, the 
amplitude of the tip of the bucket or of some other con- 
venient point, is unity. The constant a is found from the 
vibration amplitude deflection curve, which depends upon 
the type of vibration. Since every element of the integral 
in the numerator is u times that in the denominator, and 
the amplitude of every element bears a constant ratio to 
Ui, the ratio of these integrals is a linear function of ui as 
shown. Thus 


^ Trjog 

2<xwui^8 


(189) 


To find the maximum stress at some point P of section 
modulus I where the bending moment is AI, and where h 
is the distance from the neutral axis of bending to the 
extreme fiber, (189) is used in combination with the 
relation 

= (190) 


This relation is obtained from (20) of §21, remembering 
that, when y = h, a — cr„. and that l/p= dhifdx^. The 
term {dhi)/{dx-) in (190) may then be expressed in terms 
of 7<i when the deflection curve is known, and (190) sub- 
stituted in (189) in place of uy, thus giving an expression 
for the required o-„. 
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88. Application. — ^The use of this method will be shown 
by a relatively simple application : 

Consider a turbine bucket of uniform cross-sectional 
area to be acted upon by a pulsating force of the type 
assumed, which acts in the plane of vibration as shown in 


Fig. 80, and which has a constant 
amplitude of /o per unit length of the 
bucket along its entire length. The 
bucket is assumed in this case not to 
be held at the tip by a shroud band. 



Fig. 81. — Turbine Bucket 
Section. 


but to be vibrating at resonance with the exciting force 


as a simple cantilever beam. 


Other constants are as follows: 


Length / = 5-Finches. 

I of cross-section = 0.0037 inch"*. 

Distance h to neutral axis = 0.25 inch. (See Fig. 
81.) 

Decrement 5 = 0.01. 


Find the maximum fiber stress at the root of the bucket. 
The amplitude deflection curve for this type of vibration is 
found from the general expression (91) of § 39 for flexural 
vibrations of a bar of uniform cross-section. By the use 
of the conditions (93) in this same section, the constants 
A, B, C, and D can be evaluated for the case of a fixed- 
free bar, using also the condition that u = Ui when x” = /. 
Equation (91) then becomes 

u = ui (0.50128 cosh m x — 0.36806 sinh m .v] 
—0.50128 cos m x -f 0.36806 sin in x) 

where for a fixed-free bar m = 1.875// (§39). (This is 
the equation of the amplitude curve of Fig. 80.) 

^ The constant a is found by graphical integrations using 
this deflection curve, and is approximately equal to 0.631. 
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The curvature of the bucket at its root when deflected 
according to equation (191) equals 


0^ = 1.00256«F«1 . 

. . (192) 

Thus from (190) 


= 1.00256£/wm2?4i . . 

. . (193) 

Furthermore, for a fixed-free bar (§ 38), 


"’=”“(^0 • ■ • 

. . (194) 

Substituting (189) and (194) in (193), 


1.00267r/ok 

2tn^al 5 

. . (195) 


Since in = — and a = 0.631, 

<r„ = 0.71-^' (196) 

Substituting for h, I, I, and 5 the values given at the 
beginning of the section: 

(T„ = 120,000/o 

It is thus seen that if the total pulsating force (fo X 0 has 
an amplitude of 1 pound, 

120,000 , . ,2 

(T„ = — — = 23,400 pounds per mch^ 

This is dangerous from the viewpoint of fatigue, especially 
in the presence of steam. 

When the bucket is held by a shroud the resonant force 
amplitude required to produce dangerous stresses (see 
Fig. 82) is considerably greater. This case can be handled 
by the same method, using the proper deflection curve, 
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derived from equation (91), and of course calculating a new 
value of a and expressing the integrals in terms of the 
deflection at some other point of the bucket 
than the tip. There are two zones of 
maximum stress, one at the section A, and 
the other at B, Fig. 82, for which the 
cuiA’ature may be found in terms of the 
same point on the bucket chosen for the 
evaluation of a. 

Note from relation (195) that cr„ de- 
creases directly as the decrement constant 
6 is increased. 

This type of analysis tells us that 
short steam turbine buckets which occa- 
sionally fail require a pulsating force of 

^ X’ X\J>, 04.. J.U1.UUIO 

considerable magnitude to break them. Bucket Vibration 
even when the bucket vibrates exactly in withTipSupported. 
resonance with the force. This force is 
ver}”- likely to be the impulsive force on the bucket from 
passing equally spaced nozzle jets at high speed. 

89. Dynamic Vibration Absorber. — If an objectionable 
\dbration has great constancy of frequency, a powerful 
periodic force opposing this vibration may be obtained 
from a separate vibrating system which resonates to it. 

Such a de\dce is called a dynamic vibration absorber, Emd 
has been successfully used in some applications. It must 
be used vdth care, however, as a small change of frequency 
may cause it to aid instead of to oppose the objectionable 
vibration. (See reference 23.) 
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